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SYNOPSIS 


NUMERICAL SIMULATION OP TIDES AND STORM SURGES 
IN THE BAY OF BENGAL 

Vw B, K* KUMAR 

Department of Mathematics 
Indian Institute of Technology^ Kanpur, INDIA 

September, 1983 


The Bay of Bengal, v^ich accounts for nearly of the 
tropical cyclones of the Globe, is well known for the devas- 
tating storm surges* The recent exanples are the surges 
of 1971 (Frank and Hu«saln, 1971) and of 1977 (Johns et al», 
1981), With the funnel shape of its coast-line and the 
vast extent of the continental-shelf at its head, the Bay of 
Bengal is a potential region for huge tidal ranges and 
inundating surges. In spite of this fact, there are relati- 
vely few studies on the Object of tides and storm surges 
for this region as ccst^ared to the studies for the other 
regions such as the North Sea,. Motivated by this fact, an 
atten^t has been made to investigate the dynamical aspects ' 
of the tides and stoarni surges, in a detailed form, for the 
region of the Bay of Bengal through numerical si^lat^n. 

The present dissertation is presented in five chapters. 



Chapter O contains five sections. Sections 1 and 2 

I 

present general introduction and a brief historical account 
on the study of stom surges in general and for the Bay of 
Bengal region in particular. Formulation of hydrodynamical 
equations recjuired for the two-dimensional surges and tides 
is presented in section 3 and the corresponding formulation 
required for the three-dimensional surge is presented in 
section 4. But the main hydrodynamical formulation required 
for the simulation of M^-tide in the Bay of Bengal is presented 
in Chapter 1. The layout and conrtibution of the present 
dissertation are given in section 5. 

In Chapter 1, the tide has been simulated using the 
well-loiown Proudman's Tidal Theorem (Proudman 1925). Owing 
to its practical limitations (Doodson 1958, Hendershott and 
Munk 1970, Cartwright 1978, 1980) only one major application 
of the theorem has appeared in the literature - after a long 
gap of 28 years (Fairbaim, 1954). Fairbaim applied the 
tidal theorem to detemlne the K^-constltuent along the 
equator in the Indian Ocean. His selection of the constituent 
K^, owing to its frequency coinciding with the frequency of 
the earth's rotation, made the mathematical treatment of tbe 
problem simple* In the present study, however, the 
constituent is considered, as it is the major tide-producing 
constituent. Because of this, and also due to the dimensions 
of the Bay of Bengal which has a relatively lai^e meridional 



extent as compared to the extent along the parallels, the 
auxiliary functions have shown a relatively complex behaviour, 
requiring a good deal of effort in the computation of the 
M^-tide distribution. The present formulation has allowed a 
d^th variation in the meridional direction in contrast 
to the Pairbaim's selection of constant depth throughout. 

The distribution of the constituent is determined for 

lO i 

each of the spaced parallels apart) and a co-tidal map is 
presented for the region of the Bay of Bengal from 10°N to 
20®N, which excludes the shelf region at its head,. ; 

Surges, tides and their interaction are studied in I 

' I 

Chapter 2 using the vertically averaged model given by | 

i 

Plather and Heaps (1975), In the first section the numerical i 
scheme is described and surges without the tidal forcing are 
simulated, Ihe selected features include the incorporation j 

of non-linear advection terms, fairly realistic bathymetry, 1 

[ 

storms with asyronetric wind field, three types of trateks , i 

striking the coast normally and obliquely - all having typical 
features representative of the Bay of Bengal region, The 
model has shown different response features to different 
meteorological forcing features, particularly the tracks. 

Two unexpected features observed from the results of the 
numerical experiments are the swift, local oscillation, fall— 
rise-fall, of the elevation during houirs of land-fall along 
the track I II and the maximum surge peak to the left of the 



storm track 1. Both these factors have been attributed to 
the possible combination of the coastal geometry, bottom 
and the moving wind-stress field. 

In the second section a tidal regime has been established 
by forcing a sinusoidal tidal oscillation of 1 meter amplitude, 
in phase, along the open boundary. The established tidal 
regime, after 28 cycles of integration, is fairly in close 
agreement with the observed tides and is intended for the 
simulation of tide-surge interaction in section 3,. 

In section 3, tide-surge interaction, with an imposed 
initial tidal distribution of the established tidal regime 
of the previous section, is studied. Pour phases corres- 
ponding to maximum- falling zero-minimum- rising zero along 
the open boundary, of the tidal cycle have been chosen for 
this purpose and the obtained results are discussed. 

In Chapter 3 a spectral method given by Heaps (1975) 
is chosen to simulate the three-dimensional features of the 
storm surge, A selected number of numerical experiments 
have been performed with the kind of meteorological forcing 
as has been taken in Chapter 2, Regarding elevations, to a 
large extent, the major features of the two-dimensional surge 
of the previous chapter are reproduced. However, some 
interesting features of the residual circulation which can 
not be realized in a two-dimensional model are presented. 



The present dissertation concludes with Chapter A, 
giving a comprehensive review of the results obtained and 
the models used in the earlier chapters* A future program 
to improve upon the present work has been suggested. 



CHAPTER 0 


INTRODUCTION 


Roll on, thou deep and dark-blue Ocean, roll ! 
Ten thousand fleets sweep over thee in vain; 
Man marks the earth with ruin, his control 
Stops with the shore; 


-BYRON, Childe Harold, (iv). 


Perhaps ocean is one of the most fascinating features 
of the Vature that rrenkind has experienced. His experience 
has progressed through mere curiosity and helplessness 
towards understanding and harnessing the ocean phencxnena 
for the benefit of mankind and forewarning the calamities 
to save life and property. 

Among the destructive phenomena of the sea, storm 
surge is one of the most devastating* Bie fury of the 
hurricane wind causes devastating surges inundating the 
low-lying coastal areas, toy of Bengal, which accounts 
for nearly 15?^ of the world crop of tropical cyclones, is 
an excellent exaicple of this phenomenon. The following 
are the known events of the destruction caused by st«5m 


surges . 



Table 0,1 


Year 

Place 

Destruction caused by surge 

1717 

Months o£^ Hoogly River, 
near CaUcnitta 

Forty- foot Strom wave killed 
3 00,CX)0 

1364 

Months o£' Hoogly River 
near Calcutta 

50,000 were drowned 

1376 

Baclcerr Qi age, India 

Storm tide killed more 
than 100,000 

19 60 

Baia scrs; India 

About 5,000 people were 
killed 

1970 

Chit taLgaa-g, Bangladesh 

More than 300,000 were 
killed 

1977 

Di’wl Poia t, India 

More than 9 ,000 were (cilledl 


Tbe destiac-tion is so appalling and the facts are so 
appeOLtg -tlnat: oine feels the immediate need of studying and 
forewarti-ng •th.e devastation for the region of Bay of Bengal, 
dm this encSearoc recently some progress has been made in the 
stuc3y od s to 30 a surges, with the help of numerical models, in 
India aaidl Bangladesh. We describe some 6f these works in 
the naxrt sec-tioa. In the following paragraphs we briefly 
describe thae imaptortant factors involved in the surge dynamics. 

De fin it ion ol£ Surge * 

Tliere oar be several ways of defining the storm surge 
pheraonneina . Hear® we state die definition given by Heaps (1967) 





It is defined as^a raising or a lowering of sea level by 
wind and by changes in atmospheric pressure over the sea 
associated with a storm' , 

The forcing due to stmom is the principal cause for 
storm surges. There are other factors as well# which 
influence the surge dynamics significantly. Thus, the main 
factors involved in the surge generation are the intensity# 
size# speed (vector motion) of the cyclone} the track of 
the stormj the bottom topography and the coastal gecanetry 
of the sea basin# a.nd the phase of the astronomical tide. 
These are well documented by Dunn and Miller (1964). Other 
factors include the presence of river mouths and Islands. 
Local effects such as short -period waves# swell and resonance 
are also important factors in surge amplification. As an 
introduction of the subject we briefly describe these aspects 
in the following paragraphs. 

Forcing due to Storm 

Cyclone impinging on a coast provides two forcing 
components# namely# the winds and the pressure drop. Wind 
stresses (due to winds) acting over the sea surface cause a 
gradual piling up of the sea surface towairds the coast# 
while the Inverted barometer effect (due to pressure drop 
with respect normal sea level pressure) can cause an almost 
instantaneous elevation ©f the sea surface - usually 



1 cm per 1 mb (millibar) drop of pressure* However, the 
elevation due to inverted barometer effect seldom exceeds 
one meter. The size and intensity of the storm have a 
direct impact on the extent and magnitude of the piling up 
action by the associated winds. The vector motion of the 
storm, on the other hand, determines the duration of the 
storm forcing. 

Storm Track and Coastal Gecanetrv 

Storm track is an impoirtant factor in the prediction 
of storm surges. In the normal circiorastances the peak 
surge occurs to the right of the direction of the storm 
movment, because of the stronger winds there* However 
occasionally the surges also occur to the left of the storm 
track. There can be several possibilities leading to this 
type of unusual behaviour. One of th^ is presence of an 
enclosed type of coast facing the off-shore winds. The 
angle of attack of the storm, which is determined by the 
track is an essential parameter in the storm surge dynamics. 
Coastal geometry affects the surge propagation in many ways. 

A converging bay causes surge amplification through funneling 
effect. Estuaries and river mouths ^re such regions vhere 
the surge amplification leads to devastation in the adjoining 
low-lying areas. 



Bottom Topography 


Bottom topography is also a governing factor in the 

surge dynamics. In 2-D models/ the hydrodynamical eguations 

Inversely 

imply that the surge height is ^ proportional to the depth of 
water. This implies amplification of surges in shallow 
regions. Also, in shallow regions the (surface) flow 
follows the bottom contours and the geometry of bottom 
contours plays an important role in the redistribution of 
the accumulated water by the wind stresses. Consequently 
the surge prediction may be affected regarding the time of 
arrival/ the magnitude and location of the peak surges. 

The Astronomical Tide 

The astronomical tide, unlike the storm surge/ is a 
regular phenomencn associated with the periodic rise and fall 
of the ocean surface due to the presence of relatively 
moving celestial bodies. In certain coastal areas the 
amplitude of the tide is substantial. The phase of the 
astronomical tide determines the suirge heights near the 
coastal areas. If the peak surge near any coastal location 
coincides with the high tide at the location# the surge can 
be more devastating. 

Rivers and Islands 

River discharge tends bo raise the mean sea level near 
the mouths and the adjoining estuary regions, Also/ the 



rainfall associated with the striking storm may aggravate 
the situation - leading to excessive river discharge. On 
the other hand the presence of Islands may produce local 
changes in the surge propagation. 

It is very interesting to note that the Bay of Bengal 
has all the above mentioned features and provides a dwelling 
place for the storm surges. The extensive shelf area towards 
the north and the bottom topography with the presence of a 
large canyon^ the enclosed type of coasts, the mouths of 
Gangetic river system and the Islands are the remarkable 
features. 

The evolution process of surge can be divided into 
4 different stages. They are (i) the swell (ii) the fore- 
runner (iii) the surge and (iv) the resurgence. Swell 
corresponds to the rough sea conditions in the core— wind 
region of a cyclone and on the other hand 'fore-runner' 
is a slowly accumulated elevation ahead of the core-wind 
region. As the storm approaches the landfall point, the 
elevation takes place at a ripid rate and this stage is 
known as the 'surge'. When the storm crosses the coast, 
it weakens and the direct action of the storm on the sea 
surface is absent. However, in the process of returning 
to normal sea-level state, scame kind of seiches take place 
and this stage is known as ■'resurgence'. 



In the following section we present a brief historical 
account on the sub j ect of surges • Regairding tides we 

present the historical account in chapter 1* 

0*2 A Brief Historical Account 

The tulk of the literature on storm surges has appeared 
after the year 1950, However# the contributions prior to 
the World War II are quite significant and paved the way for 
greater understanding of the surge dynamics. The theoretical 
aspects of tides and waves were given in an article by Airy 

! 

which marks a beginning towards the study of storm surges, A 
definite demarcation between the early publications and the \ 
more recent publications on the subject is given by | 

Bretschneider (1967), There is a limited literature prior 
to the year 1900, Between the years 1900 and 1930 some 
important contributions towards the theory have come up from 
the authors like Abbot# Taylor# Cline# Pettersson# Jeffreys# j 

i 

Marmer# Proudman# Doodson# Yamagutl and perhaps a few others. 

The literature upto 1930 accounts for only IC^ of the bulk 
of the literature until the present time. From 1930 to 1945 
only a limited number of publications have appeared# perhaps 
due to second World War, The notable authors during this 
period are Montgomery# Rossby# Arakawa and Yoshitake, An 
extensive bibliography on storm surges upto 1954 is given by 
Kramer (1955), 



1950s 

Perhaps more than 80% of the bulk of the literature 
has appeared since 1950. With the advent of large^ high- 
speed computers the literature began to expand at an increasing 
rate. Three particular regions in which extensive studies 
have been made are the coastal regions of United States, the 
North Sea and Japan, The works of Goldsbrough (1952), 

Proudman (1953, 1955), Lauwerier (1955), Crease (1956), 

Doodson (1956), Hansen (1956), Reid (1956, 1957), Harris 
(1957, 1959), Platzman (1958), Weenink (1958), Kajiura (1959), 
Rossiter (1959) and Wilson (1959) are some of the notable 
contributions that have appeared during the decade of 1950s, 


1960s I 

An excellent article on the numerical prediction 

i: 

of storm surges is given by Welander (1961), A summary of 

i 

all the available, pertinent information on the state of the 
art of the subject on storm surge (up to 1965) is given by | 

Bretschneider et al. (1966), Ihis is followed by an extensive | 
bibliography on storm surges by Bretschneider and Pick (1966) 
and a review by Bretschneider (1967), Modeling is one of 
the dimensions acquired by the storm surge studies during 
this period. Theoretical models have been given by Heaps 
(1965) to investigate the surge on a continental shelf and 
by Heaps and Ramsbottcm (1966) to study the wind effects in 
a long narrow lake. With the development large— memory, 
fast-computing machines of later generations, n\americal models 
have been developed 'to include more complex problems, A two- 



dimensional numerical model has been developed by Heaps (1969)# 
incorporating the general sea boundaries within the basic 
grid system. A brief review of some of the numerical schemes 
can be found in this work. Many of the numerical schemes 
used by them in solving the governing equations by the 
explicit finite— difference methods. Seme Implicit schemes 
are given by Leendertse (1967). Other Important works during 
this decade are the numerical models by Jelesnianski (1965# 
1966)# Miyazaki (1965). A general descriptive article on 
surges is given by Groen and Groves (1962). 

1970s 

A further leap in the early seventies is the 
development of three-dimensional models. Some of the models 
which belong to this class are multi-layer and multi-level 
models which are discussed in a review by Liu and Leendertse 
(1978). Three-dimensional models have also been reviewed 
by Davies (1981). 

A new approach to the three-dimensional modeling for 
the study of storm surges is given by Heaps (1971). This is 
a spectral method in which the horizontal components of 
current are expanded in terms of eigen functions developed 
through the depth coordinate. Coefficients of the expansion# 
varying in the horizontal and with time are determined f ran a 
two-dimensional numerical scheme. The developing three- 
dimensional current structure# and the elevation patterns are 



determined by means of converging series involving these 
coefficients and the depth coordinate. With the develojxnent 
of this model, a number of applications are given by Heaps 
and his colleagues, A three-dimensional numerical model 
of the Irish Sea is given by Heaps (1973, 1974) and Heaps 
and Jones (1975), As a further development, a formulation 
of the three-dimensional model, to include nonlinear terms, 
is given by Heaps (1976), The spectral model has also been 
used to determine the density currents in the Irish Sea by 
Heaps and Jones (1977), and to determine tides and storm 
surges with vertical eddy viscosity prescribed in two layers 
by Heaps and Jones (1981), During this period three- 
dimensional numerical models have also been developed by 
Forristall (1974), Durance (1976), Davies (1977(a), (b)) 
and Nihoul (1977), and others. Some important two-ddmensional 
models have also been developed during the decade. Among 
them are the works by Banks (1974), Prandle (1975), Prandle 
and Wolf (1978) regarding tide— surge interaction} Jelesnianski 
(1970, 1972, 1976), Reid et al, (1977) , In 1978 World 
Meteorological Organization has published a repoirt on 
tropical storm surge prediction which gives the present day 
art in surge prediction. A comparative assessment of the 
numerical models for the North Sea has been given by Nihoul 
and Ronday (1976), Sesides these works there are several 
other important works vhich have made significant contributions 
to the understanding of the storm surge dynamics. 



study of Stom Surges for the Region of the Bay of Bengal 


Storm surge studies for the region of the Bay of 
Bengal are quite recent as compared to the studies for the 
coasts of United States and the Horth Sea* The initiation 
has come from the articles by Frank and Hussain (1971)^ 

Das (1972) and Fleirl and Robinson (1972)* Das et al. (1974) 
have investigated the surge effects due to the storm that 
struck Chittagong (Bangladesh) in 1971. A two-dimensional 
formulation is given by Johns (1981) regarding the simulation 
of surges in the Bay of Bengal* This formulation has been 
used in a subsequent study by Johns and Ali (1981) for the 
coasts of Bangladesh. A multi-level model with a-coordinate 
approach is also given by Johns (1981)* Recently Das (1981) 
has studied the tanporal response of wind stresses and has 
shown that in the initial stages of the surge the diveirgence 
of the wind stress is more Important than the curl of the wind 
stress# while in the later stages the curl becomes more 
important* Ihis response has been related to the storm 
speed (vector motion) and it has been shown that the response 
is inversely propoirtional to the square- root of the basin depth 
An application of this model has been given by Sinha (1931) 
for the Gulf of Thailand* For the east coast of India# surges 
are studied by Johns et al* (1981;) • Besides these works a 
few other works have appeared regarding some individual aspects 
of the surge for the region of Bay of Bengal. Of these# the 




studies by All (1980), Dube et al. (1981), Johns et al. 
(1982) are notable* 

In the next section we proceed to describe the 
hydrodynamical formulation for the study of surges and 
tides *. 


0*3 Hydrodynamical Equations and formulation of Two- 
Dimensional Surge 

We staart with a more general form of hydrodynamical 
equations of Reynolds type to represent an ocean on the 
rotating earth* We consider the right handed system of 
coordinates ‘ With the z-axis upwards and the xy plane 
coinciding with the mean sea level, xj-axis directed towards 
the east and y-axis directed towards the north* 


With these considerations, the hydrodynamical equations 
under the action of a tidal potential q' are 


at ^ ^ ax ^ ^ 3y ^ ^ 


f V = - i aE » + 

^ p ax ax ^ 


x i. r T X T X t 1 

P ^ ax XX 3y yx az ^zx J 


(0.3.1) 


at ax ay az 


P ay ay 


Ir X ^ X i- t X 1. T 
P ^ ax xy ay yy 32 zy 


(0*3*2) 



g 


3w 
3 t 


+ u 


3 w 
3x 


+ V 


dJU 

3Y 


+ w 


32 ■ 32 


+ i.ri-.T +^T +^T 1 

P ^ 3X X2 ^ ay yz ^ 3z zz -« 


(0.3,3) 


and the equation of continuity is given by 


an + az +2W ^ 0 , 

3X ^ ay ^ 3z ^ * 


(0..3.4) 


Based on the scales of motion we have to simplify the above 
equations to allow a solution of them for a given set of 
initial and boundary conditions. 

In the study of tides and surges it is usual practice 
to approximate the equations (0.3.1) to (0*3.4) by means of 
long ocean waves for which it is assumed that the free 
surface elevation is small and the horizontal scales of 
motion are large as compared to the depth. With the help 
of dimensional analysis (Welander 1961) one can neglect 
vertical accelerations and vertical convection terms# 
justifying the hydrostatic approximation in the vertical. 
Under these assumptions equation (0*3,3) becomes 


or 


3z 


gp 


P » Pg - gp (z - ^ ). 


(0,3.5) 



and 


We also ass-ume that the lateral stresses T it 

xx' xy 

Tyy to be insignificant for the present investigation as 
we are not considering narrow regions for the adjacent 
currents to be important (Groen and Groves 1962)4 Under 
these asstimptions^ equation (0,3*3) can be eliminated by- 
incorporating its reduced form (0.3.5) in the equations 
(0,3*1) and (0,3*2), Thus we have 
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where we have replaced the tidal potential q* by -the 
instantaneous response (called -the equilibrium elevation) 
and for convenience it is written 


ail = _ L_ 

3x 3x pg ax 

ail =iL - i- 

ay ay pg ay 


(0.3.9) 


The common type of boundary conditions associated 
wi-th the above hydrodynamic al system of equations 
(Welander 1961) are 



u = V = w = 0 


T = F, / T = G, 
zx zy b 


at the bottcxn z = -h 


and 


(0.3.10) 


dl Zl 
^ ^ dt ^ 3t 
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I at the free surface (0.3.11) 


In the actual problem solving, the hydrodynamical j 

equations are subjected to appropriate initial and lateral f 
boundary conditions^ in addition to the boundary conditions 
(0.3.10) and (0.3.11). 

For convenience in the later stages of our study we 
depart from the present system of equations in two different 
directions. In the first instance we derive a system i 

of equations that is amenable to the study of two-dimensional 
surge. In the second instance we derive a systan of 
equations leading to the study of three-dimensional surge. 

I 

Formulation of the Problem for the Study of 
T wo-Dimensional Surge 

In an attempt to simplify further the system (0.3.6) 
to (0.3.11) it is usual to introduce the vertical averaging 
defined by 

I ^(x^yjt) ^ 

(x,yrft) = 7 ^ / u(x.y, z,t)dz . (0.3.12) 

^ -h(x#y) 


With this vertical averaging and the use of boundary 
conditions (0,3,10) and (0,3 ,11) the systapa (0.3*€) to 
( .0 „ 3 *8 ) becomes 

3 

+ ~ { (S +h) (uv - u V ) }] - fv 

= -g ^ (.l-V) + (P3-Pb) ( 0 , 3 . 13 ) 

tIkC h - Sv)) 

+ ~ £(5+h)(v^ - V )}] + :^ 

= - g ^ ( 5 -S') + 

|| + f_ j(5+h)5) 

where 

F = T 1 F, 

s zxl Z= £, , b 

G = X I G, 

s zyl z= <, , b 


+ ^ {(S+h) V} =0 (o.sas) 

dY 

” ^zxl z= —h 
“ ^zyl z= ~h 



(0.3.16) 



Though the above equations are vertically averaged^ 
some form of vertical structure still prevails in the sense 
that the square of the mean (product of the mean) and the 
mean of the square (mean of the product) for the variables 
u and V need not be the same* In the absence of 
extremely shallow water this vertical structure can be 
neglected with the assumption that u = u^ , v = v^ and 
uv = u V. Also^ in the expressions for and 2;^ we 

do not consider the equilibrium tide and instead we enforce 
the tide as an input from the adjacent oceans through the 
lateral open boundaries* Thus, in the present case we 
have 


all _ 1 

3x - pg 8x 


all 1 

ay pg 9y 




(0.3.17) 


Taking the above arguments into consideration the 
equations (0.3.13) to (0.3.15) with the over-bars removed 
can be written as 


3u . .. au . au ..a? 1 9^a . ^b 

+ - p^TST 


**■ ^ + V ^ + fu .= -g i 


at 


ax 


dy 


By ~ p ay 'P(t+h) 


It + fe ^ ^ {(S+h)v} = 0 . 


(0.3.18) 

(0.3.19) 

(0.3.20) 



The closure of the above system requires the characteri- 
zation of the bottom-stress components in terms of 

the mean flow variables u^v* Depending upon the t]^e of 
flow in the model we can have relations between the bottom 
stresses and the mean flow quantities, A good account on 
the subject can be found in Welander (1961), In the present 
work we adopt the widely used quadratic law for bottom stress 
characterization. Thus 

= kPu (u^ + = 

= kPv (u2 + v2)V2 ^ 

where X = kp (u^ + 


K u 

Xv 


(0.3.21) 


1 

In addition to the parameterization of the bottom 
stresses the problem requires the specification of stresses 
near the free surface as an input fix>m the meteorological 
forcing. This is provided in terms of wind stresses near 
the free surface. Though not directly concerned with the 
ocean dynamics (the present problai:i) we have a two— stage 
oroblOT in the specification of the wind stresses. In the 
first stage one should have a good observed values of the 
surface winds and in the second stage one should choose 
i good parameterization of the wind stress^ (at the sea 
surface) in terms of the observed wind velocities. Here too. 



we can have a variety of relations in practice depending 
on the type of the approximations made for the surface winds 
from the observed winds^ usually 10 meters above the sea- 
level* A good account of these aspects can be found in the 
review by Lick (1975)* As in the case of bottom stresses 
we adopt the commonly used quadratic law for surface 
stresses given by 

(0.3.22) 

where w refers to the wind condition. 

Other boundaries to be included in the model are the 
lateral boundaries of closed and open type. For obtaining 
a solution of the equations we have to prescribe conditions 
across these boundaries. In the present formulation we 
consider the fixed wall-type closed boundaries across which 
the water transport is set to zero. In real situations 
incorporation of moving coastal boundaries would be proper. 
But in the absence of reliable data regarding the coastal 
topography one can consider the wall- type coastal boundaries 
in which case, the obtained elevation magnitudes, for obvious 
reason, tend to be higher. Regarding open lateral boundaries 
we can prescribe various types of the conditions. In the 
present case we prescribe the following conditions for the 
lateral boiindaries* 
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(0.3.23) 


For coastal boundaries : u.n = 0 

For open boundaries : S = 0 in the case of 

no- tide 

or S = A cos at in the 

case of tide with 
amplitude A and 
frequency a 

(0.3.24) 

To avoid the reflection of long a waves from the models’' 
interior special procedure has to be used at the model's 
open boundary. In an analytical investigation of coastal 
tides Proudman (1940) applied the admittance type of 
radiative boundary condition. The condition has also been 
considered by R&dd and Bodine (1968). Masch et al. (1971) 
converted the forcing tide at the model's open boundary into 
a tidal flow according to the relation 

^n = ^'^g ‘ (0.3.25) 

where C is an admittance coefficient taken as the 
gravity wave velocity, ^ represents the specified water 
level at the boundary (eg. observed time history of water 
levels) and ^ is the previously computed water level at 
the boundary. 

In view of the characterization of bottom stresses we 
rewrite the equations (0.3,18) to (0.3.20) , with the continuity 
equation appearing first. 



Thus 


II + ^ {(?+h)u} + ^ {(2+h)v} = O 


(0.3.26) 


at 3x 3y 



1 Xu 

p air "jirffi) + HTShT 


(0,3.27) 


2v 

3t 


+ u 


av 

3x 


+ V ^ + fu 

3y 



XV _fs 

P ay strafi pg^j,) 


(0,3.28) 


In chapter 2 we directly incorporate these equations 
for the numerical treatment along with the appropriate 
initial and boundary conditions. 


0 ,4 Formulation of the Problem for Three Dimensional Surge 

Three-dimensional models for tides and storm 
surges are physically more relistic as compared to the two- 
dimensional models. Vertically averaged models are simple 
for numerical treatment and a good number of e 3 q)erimentatlons 
can be performed at a reasonable effort but forgo some of 
the physical processes required for better undeirs tending 
of the surge mechanism. For example^ in a two-dimensional 
(vertically averaged) model the characterization of the 
bottom friction^ for conceiving the 'closure"^ # is done 
dn terms of the averaged velocity components# instead of 
actual velocity components near the bottcm. 


Ihe three- 



dimensional model can closely represent the current structure 
whereas it is not possible in the case of two-dimensional 
models. However# the three-dimensional models are much 
more tedious and require enormous computer budget as well as 
human effort as compared to their two-dimensional counter- 
parts* Perhaps a better way is to have a compromise 
between the two-dimensional and tlio three-dimensional models 
by splitting the complex probltan in terms of tasks moulded 
into either of these models# complementing each other# in 
achieving the required objective. 

Some notable three-dimensional models have been 
developed in recent years. Heaps (1972# 1973# 1974# 1976) 
has developed a spectral method for the numerical solution 
of the three-dimensional hydrodynamical equations for tides 
and storm surges* In this method# the horizontal components 
of current are expanded through the vertical frcan the sea- 
surface to the sea-bed in terms of a set of eigen functions# 
taking into consideration# the dynamical boundary conditions 
at the free surface and at the sea-bed. Coefficients in 
the expansion# varying in the horizontal and with time are 
determined following which the current structure and the _ _ _ 
elevation pattern can be determined. Application of this 
method is given by Heaps and Jones (1975# 1977). In all 
these models# however# a constant eddy viscosity throu^ 
the vertical is assumed with a slip condition at the bed 



representing^ albeit crudely# the special circiomstances 
hold in the bottom layer* Nihoul (1977) has developed a 
three-dimensional spectral model# some vjhat similar to 
that of Heaps* model but incorporating a variable eddy 
viscosity in the vertical* In the later models, however. Heaps ( 
Heaps and Jones (1981) have considered linearly varying 

vertical eddy viscosity in their two-layer model. In this 
model they have resolved a bottom boundary layer in vdiich a 
linearly varying eddy viscosity in the vertical ( logairithmic 
velocity profile) is taken with the assumption of no-slip 
at the sea-bed* In the upper layer a constant vertical 
eddy viscosity is assumed making the distribution contimious 
at the interface. Johns (1977) has developed a three- 
dimensional multi-layer model with a-coordinate approach* 

An application of this is given by Johns et al* (1982)* 

Other notable models are given by Davies (1977 a,b), 

Porristall (1974), and Durance (1976)* Despite these works, 
the developanent of three-dimensional nxunerical models of 
water movements in the continental- shelf sea is in its early 

!■ 

stages and much has to be done in developing the models with j 
vertical density stratification so that thermocline movements 
and the associated currents may be taken into account and 
predicted (Heaps 1981a) * In the following paragraphs we 
describe the spectral model given by Heaps (1976)* j 



Rewriting equations (0.2.6) to (0.2*8) in a convenient 
form# we have 


aL!i+u2Lii4.v2Ji + w2^-fv = 

at ^ ax ^ 3y ^ 3z 


-g (S- SO +1^ (N ||) 
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/ udz + 

-h 



vdz = 0 


(0.4,2) 

(0.4.3) 


in which the continuity equation is retained in the 
vertically integrated form. In writing the above equations 
we have assximed the eddy- viscosity model with N as the 
coefficient of vertical eddy viscosity. These equations are 
to be solved for u^v and t subject to 

u=v = w = 0 at the bottom z = -h(x/y) (0,4.4) 

and 

W = —sr =s r-r + U ^ + V II: at the free surface 
at 9t 9x 9y 

z = S (x/y/t) (0.4.5) 

and also subject to appropriate initial conditions and 
lateral boundary conditions. 

The vertical velocity ccanponent is given by 

w = ^ / udz +1“ J vdz . 

-h -h 


(0.4.6) 



By the choice of the coordinate system the free 

surface of the ocean is at z = g (x/y^t) and the bottom is 

at z = -h(x,y). If P * G denote the x^y components 

s s 

of the wind stress on the sea surface/^ then 



( 0 , 4 . 7 ) 


Correspondingly, if denote the x,y components 

of the frictional stress of the vrater on the sea bed, then 


^b = 


|^>-h 


G^ = - P (N 
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z''— h 


( 0 . 4 . 8 ) 



( 0 , 4 . 11 ) 



In addition to these conditions we consider the 
problem in a more physically realistic sense. For example, 
the non-linearities result from the occurence of ^ as 
the lower limit in the integrals of equation (0,4*3) * due 
to the presence of the convective terms in equations 
(0,4,1) and (0,4,2); because of the evaluation at 2 = 5 lin 
the equation (0,4.7); and because K is a function of the 
motion in equation (0,4,10) , We include all is facts in the 
formulation and any modification would be mentioned instantly 
when such an occasion arises. 

Trans f omiat ion 


In the present method (Heaps 1976) we introduce a 
function (p(x,y, z,t) defined within the depth range 
$ > z > -h for each x,y/t. Thus 


(p = (p(x,y,z,t) “h < z < t (0,4,12) 

where 

I = I (x,y,t), h = h(x,y). (0.4.13) 


Multiplying the equations (0,4.1) and (0,4,2) by <p(x,y# 2 ,t), 
integrating from -h to t and then dividing by $ + h, 
yields 
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where 
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The terms and arise from the following operation 

°° H ^“'5 li • 




(0.4.21) 




(0.4.22) 


The terms and arise from the convective accelerations . 

The vertical structure is retained and is represented by means 



of a coipplete orthonormal set of functions appropriate to 
the present problem. 

We observe that the equations (0.4.14) and (0.4,15) 
have some comparable similarities with the equations (0,3,5) 
and (0,3.6) respectively. In this perspective, it is natural 
to loo]c forward for the characterization of the integrals 

" h .y h 

to make the equations (0.4.14) and (0.4.15) into the familiar 
solvable form of the two-dimensional model. 


In this endeavor let us consider the integral 
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Similarly for the integral J ^ rrr (N |--) dz # we have 
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Thus if we require 
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(where X is intended to play a similar role as in the 
two-dimensional problem) we shall have to encounter# for each 
(x#y#.t)# the eigen value problem 

(0.4.27) 


(0.4.28) 

a(p 

and ~ N ^ at z = -h . (0.4.29) 

The above eigen value problem is the familiar 

3 3 

problem of Strum-Liouvelle type. The operator ~ (N - 5 — I 

o Z oZ 

for the above eigen value problem is self-adjoint and 

consequently the eigen functions o # corresponding to the 

• T 

unboundedly increasing eigen values form a ccmplete 

orthonormal. set# with the normalizing factor 
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For each (x^y^t), let the ascending eigen values 

and the corresponding eigen functions derived from the 

equations (0,4.27) to (0,4.29) be denoted by and (p ^ 

(r = respectively. Thus for each pair 

( X , (p) we have from (0.4.14) and (0.4,15) 
r r 
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where u^/ v^z a^ etc, have the previously defined meaning 
with the subscript r. Solving for u^ v^ (r = Iz^zSz...#®®) 
we can write 


OO ^ CO 

u = E u cp ^ V = S V (0.4.33) 

r=l^^^ r=l^ 


Substituting (0,4,33) in (0,4,3) we get 

^ i f ^ c(h+2)a^'j-^ =° 


(0.4.34) 



Similarly (0,4*33) in (0,4.4) yields 


OO 

w= ^ t ~ { (h+S)qj_ u^} + {(h+2)q^^j. v^} ] (0.4.35) 

where 

z • 

^r ~ h+Y ^ (0,4.36) 

The depth mean values for u and v are given by 
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(0.4.37) 


In the equations (0,4.31) and (0,4.32) the terms j~ / 

R , S ^ S take the following form 
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where 
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A particular Formulationt 


For mathematical simplicity we assume the eddy viscosity 
N to be uniform with respect to the depth coordinate z* Thus 
if we use the normalized form of depth coordinate/ given by 


= ^ZJ.. 


(0.4.49) 


the eigen value problem (0,4.27) to (0.4,29) can be written as 
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Conditions at g = 0 lead to the solution 


<p = cos (ag) 


where a = ^ 


^ (h +&) 


(0,4.53) 


Condition at = -1 Implies 


a tan a - c 
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where c = (h + t) 


(0.4.54) 
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Let a^/ r = 1/2/3/.,, be the positive roots of the 
equation (0.4.54) in the ascending order. 


Then the 



corresponding Strum-Liouvelle problem leads to 
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In practice# the series for u and v in (0*4»33 ) have 

to be terminated at some point# say# at the terra* 

Similarly all the siammations are to be truncated after the 

first M terms. In this particular formulation all the 

terms vanish. Making 

r#j#m r#j#m ' r#j#m * r#j#m 

use of the equations (0.4,56) to (0,4,60) the expressions 

of (0,4,38) to (0,4,41) reduce to 
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Using the equations (0.4.56) to (0.4.69) and with the 
properly specified boundary and initial conditions we carry out 
the canputations regarding the study of three-dimensional 
surge. This will be presented in chapter 3.. The vertical 
velocity canponent can be determined from (0.4.35) and the 
depth-mean values for u and v can be obtained from (0.4.37) • 



0.5 


Layout and contribution of the present worlc 


In the following paragraphs we give the layout and the 
contribution of the present work* 

Chapter 1 deals with the study of determining the 
distribution of M^-tide in the region of the Bay of Bengal 
by means of the well known Tidal Theorem given by Proudman 
(1925)*. The Bay depth is allowed to vary in the latitudinal 
coordinate* Amplitude and phase distributions are determined 

i O 

on each of the equally spaced parallels apart) from 10®N 
to 18°N, A co-tidal chart is drawn and is found to be in 
reasonable agreement with that given by Pekeris and Accad 
(1969) . 

In chapter 2 we have taken up the study of two- 
dimensional (vertically averaged) surge# tide and their 
interaction for the region of the Bay of Bengal* A shelf - 
area at the head Bay of Bengal is considered for the numerical 
experiments with the tidal and meteorological forcings* 
Asyiranetric storm forcing is considered to study the difference 
in the surge behaviour as compared to the surge behaviour 
due to the syrranetrlc storm (from v^ich the asymmetry is 
derived) , Other features included are the geos trophic terms, 
nonlinear advective terms, the atmospheric pressure gradimit, 
quadratic formulation for the surface and bottom stresses and 
the fixed \i^ll— type coastal boundaries* In all the numerical 



e 3 <perlments a severe cyclonic storm typical in strength and 
character for the region of the Bay of Bengal is considered. 
Three types of tracks are chosen for the ntunerical experiments 
making different angles of incidence with the coast line. 

Track I is chosen to represent the nearly normal incidence and 
track II is chosen to represent the nearly parallel incidence 
at the point of landfall near Chittagong. Track III is 
chosen to represent the normal incidence at the point of land 
fall at a location to the right of Hoogly estuary. Frcsn the 
numerical experiments with the storm forcing along tracks I 
and II (with the same point of landfall ), it is shown that 
the location of maximum peak surge occurs to the right in the 
case of track II and to the left in the case of track I. It 
is also observed that the asynanetry of the wind field (with 
the stronger winds to the right of the track) can produce a 
maximum peak surge in excess of 2 raeteres as coitpared to the 
predicted surge with the symmetric storm forcing. In the 
case of track III some local*, short lived oscillations are 
obtained during the landfall hours. Tidal regime has been 
established with a sinusoidal forcing of 1 meter anplitude 
( in phase) along the open boundary of the model area. The 
obtained tidal anplitude (2.1 meters near Chittagong) as 
well as the phase is shown to be in close agreement with the 
observed values. Tide-surge interaction has be^ studied with 
the forcing of the established tidal regime and the asymmetric 



storm forcing* Uie interaction produced by the present model 
is shown to be weak and it is felt that the further refinement 
of the mesh in the shallow areas and the inclusion of further 
river-shallow systems are necessary for the proper represen- 
tation of the non-linear interaction of tide and surge* Owing 
to the incorporation of fixed wall-type coastal boundaries, 
the obtained values of elevation are found to be higher in 
magnitude (upto 20 % ) as compared to the available observations, 
particularly near the concave comer points* However, the 
qualitative features are in reasonable agreement with the 
available observations* 

In chapter 3 a spectral three-dimensional model (Heaps 
1975) has been considered for the simulation of the three- 
dimensional surge* Same features regarding the storm forcing 
as taken in chapter 2 are considered for the numerical 
experiments* Three numerical ejqperiments, with the asyiranetric 
storm forcing along tracks I, II and III have been carried out*. 
The obtained elevation patterns of the spectral model have 
shown marked similarity with the corresponding elevation 
patterns of the two-dimensional model* The elevations obtained 
in the spectral model have shown slightly higher magnitudes 
(upto 10 56 ) as compared to the corresponding elevations of 
the two-dimensional model* The local oscillation near the 
point of land fall for the track III numerical experiment in 
both the models are quite similar* A better eaqjlanatlon of 



this behaviour has been given with the help of the spectral 
model* Also, it is verified that this behaviour is due to 
a very particular combination of the storm, coast line and 
the bottom contours. For this purpose the two-dimensional 
model has been run with the tracks passing through the same 
point of landfall (of track III) but with different 
(obliquely opposite) angles of incidence* No such local, 
swift oscillations have been observed from these runs. 

Regarding the current patterns, few interesting results 
are obtained. Daring the intensification period (the first 
ten hours after the initialization) of the storm to the 
mature stage, the surface currents (just at the free surface) 
in some regions have attained the magnitudes (over 6 meters/ 
sec !?) in excess of the so far observed values upto 3% of the 
speeds of the prevailing surface winds. In the absence of 
data regarding the currents in the severe storm conditions a 
conparison could not be made. Even for moderate wind 
conditions the current meter readings, to the best of the 
author’' s knowledge, are usually available at a level 4 meters 
(or more) below the free surface, EjKiepting these local 
excessive magnitudes, for which no field observations are 
available, the over-all current structure is in reasonable 
agreCTient with ■the normal observed values, O-ther interesting 
fea-ture is -the horizontal circulation of the water mass 
across ■the open boundary owing ■bo idle wide extent of •the open 



boundary^ and the smaller depth as compared to the 'Ekman 
Depth' for the region under study (see Appendix at the end 
of chapter 4) • The inflow across the open boundary near the 
central portion is compensated by the outflow through the 
side ways across the open boundary and is in consistence 
with the horizontal circulation pattern described by 
Welander (1961) Ihis region of inflow is observed to 
occur to the south-west of the storm centre in a region 
outside the maximum wind core. 

In chapter 4 , an over view of the present study has 
been given by making comparisons of the models used and the 
results obtained. Vast scope for the future work for the 
region of the Bay of Bengal has been suggested. 



NOMENCLATURE (Chapter O) 
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(Symbols not explained in Chapter 0) 

the velocity components taken in the directions 
of increasing X/y,z respectively 
the pressure 

the density of sea water, assumed constant 
(1025 kg/Mtr^) 

elevation of the sea~surface above the undisturbed 
sea- level i 

the depth of the undisturbed sea i 

the Coriolis parameter | 

time 

the acceleration due to gravity 
the sea-level atmospheric pressure 

the horizontal forcing stress components (due to wind) 
at the sea surface 

the horizontal stress components near the sea-bed 
coefficient of wind drag 
the density of air 

the horizontal components of the v/ind velocity at 
the sea surface 

the vector with components u,v,w 

the unit outward-normal vector 

drag coefficient (for the flow at the sea-bed) 
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( )j , ( )_ 

■'xx'^x 

Units 


refer to evaluations at z = t and 
at z = ~h respectively 
are the components of the Reynolds' 
stresses of the fluid motion 

M.K.S, 



CHAPTER 1 


M^-TIDE IN IHE BAY OP BENGAL 

1 • . Introduction 

Oceanic tides^ which we observe as periodic rise and 
fall of the ocean surface in syncronization with the motions 
of the sun and the moon^ are due to the gravitational 
attraction between the earth and the relatively moving 
celestial bodies. The subject of tidal phenomenon is an 
old one, with a vast literature, mostly Victorian (Hendershott 
and Munk 1970) . Newton was the first to give a physical 
explanation of ocean tides. Later Bernoulli, Laplace, Hough, 
Airy, Kelvin, Darwin and Poincare set up a classical theory 
of tides, the aim of which was to understand qualitatively 
and quantitatively this natural phenomenon. The discussions 
on this subject for the past three centuries have passed 
through myths, paradoxes and heated controversies (Schwiderski 
1980) . Despite its age and its vast literature, the research . 
on this subject is flourishing in new directions from time to 
time. In this respect, Doodson''s review (1958) has stimulated 
work on numerical models. At present, the stress regarding the 
development of the subject is laid on the observational aspects 
to understand the pitfalls of the hitherto developed work which 
had led to some controversies (Cartwright 1977). 



45 


The theory of oceanic tides is of special Importance 
since the observations are only available from coastal areas. 

In this regard, late Prof. Proudman (1925) had enunciated 
a theorem in tidal dynamics as an integral formulation 
involving certain auxiliary functions. This formulation 
envisages that, the tidal distribution can be determined in 
the interior of an ocean region bounded by a coast line 
(which can have openings to the other adjacent ocean basins), 
when the tidal data in the form of harmonic constants is 
available along the coast line. Proudman has provided some 
illustrations for simple land-locked basins as well as for 
basins having openings to other ocean regions. 

It was after many years, a detailed application of the 
tidal theorem has been given by Pairbaim (1954) in which 
the distribution of the semi-diurnal constituent along 

the equator in the Indian Ocean has been determined. 
Consideration of the harmonic constituent natural 

because theoretical investigations prefer in the case 

of semi-diurnal tides and in the case of diurnal tides - 

in order to avoid mathematical difficulties. Fairbaim 
in his work has assumed constant depth for the Indian ocean. 
Based on the obtained tidal distribution along the equator, 
a CO— tidal map has been presented. However, it is more 
proper to find the tidal distribution along the other parallels 
of the region for presenting a co-tidal map* In fact, this 



has been mentioned by Fairbairn himself and also by Doodson 
(1958), Prom the review articles by Doodson (1958)/ 
Hendershott and Munk (1970), Hendershott (1972^ 1973, 1977)/ 
Wunsch (1967, 1972, 1975), Cartwright (1977), Leblond and 
Mysak (1979)/ Schwiderski (1980) and from a survey of the 
recent literature available on ocean tides, it is felt that 
Fairbairn* s work still stands alone regarding the major 
application of the Proudman's tidal theorem-perhaps owing 
to its practical limitations (Cartwright 1977) • 

In the present investigation, the distribution of the 
harmonic constituent (which is the principal harmonic 

constituent for the Bay of Bengal) is determined for the 
region of Bay of Bengal. In contrast to the constituent 
K^, the consideration of the constituent recfuires a 

different development of the auxiliary functions. Depth is 
allowed to vary with the latitude coordinate. This is 
fairly reasonable in the case of Bay of Bengal which has 
narrow shelfs to the east and west. The distribution of 
is determined along the parallels of 4.atitude, spaced 

10 

apart, ranging from 10®N to 20®N. The relatively 
shallow area is excluded because it requires the inclusion 
of frictional terms in the formulation and the problem 
becomes difficult to tackle. Based on the obtained 
distribution along the various parallels running through the 
bay, a co-tidal chart is drawn. 
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1 .1 Tidal Theorem 

Hydrodynamic al Equations 

VJith the assTimption of linearized shallow-water 
equations on a rotating spherical earth, the hydrodynamical 
equations are given by (Proudman 1951) 
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and the continuity equation, is given by 


( 1 . 1 . 1 ) 


( 1 . 1 . 2 ) 


On the other hand the elevation Z =ind the current 
components u,v are also set up by the harmonic constituents. 
Thus 
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where Re stands for real part. Amplitudes and phases^ for 
instance and k ^ are called harmonic constants; 

normally they are derived from long-term records of the 
sea-level. These records are also used in the derivation 
of tidal currents. The introduction of these kinds of 
harmonic constituents simplifies the mathematical treatment 
of tidal problems by means of the hydrodynamical equations, 
Thus^ with the introduction of time dependence in the form 
^lat variables u,v and S' we can write the 

equations (1.1.1) to (1.1.3) as 


icru — fv 


.g 

a cos cp 


dt' F 

3X H 


(1.1.5) 


iav + fu 



( 1 . 1 . 6 ) 


and ia S { I— (hu) + (hv cos (p ) } = 0 (1 .1 .7) 

a cos 9 Qx 


We introduce the auxiliary functions U,V^Z which 
are analogous to the elements in a free type of motion, taut 
for the rotation of the earth (taken in reverse direction) 
and require them to satisfy the following corresponding 
equations 


icJU + fV 


.g 3Z 
a cos cp 3X 


( 1.1 . 8 ) 


iav - fU = - 3 
a 


iZ 

3<P 


(1.2.9) 
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For a region of the ocean/ Stokes’' theorem reduces to 

^ i-lsr^f h co.<P)!dS (Ul.lV) 

Replacing XLfV,v by ^ 'U / S'V, ^ 'N respectively 
in the equation (1.1,11) and using the equations (1.1.5)/ 

(1.1.6) and (1.1*10)/ we get 

/ h ds = ia // I'Z dS ~ i /; h {ld(Uti+Vv) - f(Uv-Vu) .• dS 

y 

~ I //(FU+GV)dS (1.1.12) 

y 

Similarly , replacing u^V/V by Zu# ZV/ Zv respectively 
in the equation (1.1.11) and using the equations (1.1 ,8)/ 

( 1 .1 .9 ) and ( 1 *1 . 7 ) / we get 

/ hZV ds = -ia IS Z ^ dS - I // h {ia(uU+vV) + f(uV-W)}dS 

(1.1.13) 

Subtracting equation (1.1.13) from equation (1*1 .12)/ 
we get 

/ h S'N ds - S hZV ds = ia Sf f ZdS - i jj (PU t GV)dS 
or /h(S~ r)Nds - IhZv ds = ia // fzdS - i // (KJ+GV)dS. (1.1.14) 



This is the tidal theorem which was enunciated by 
Proudman in 1925. The elevation t can be determined in 
the interior of an ocean region S, bounded by the coast 
line C, with an opening L (which we specify at a later 
stage), if the equilibriim elevation f and the external 
friction terms F,G are known over the entire area S 
and the value of 1 over the boundary. The equilibrium 
elevation can be determined from the following relation 
( P roudxian 1953.) 

S = H e3q> {2i( X + o)t - ^ 7i) } cos^cp • 

If the standard meridian is set to coincide with the 
meridian of Indian Standard Time (i.e., 82.5®E), this 
relation becomes 

f = H exp {2i ( X - I 7T) } cos^ <p (1.3 

where the value of the amplitude H for the present case 
is taken as 0.267 meters. 

The region of Bay of Bengal is considered as a region 
bounded by the coastline and open boundaries as shown in 
Pig, IfiL. We exclude the relatively shallow region above 
20®N whose study requires a separate investigation because 
of the involvement of the frictional teims F and G, For 
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an ocean region such as this, the tidal theorem (1,1 *14) 
can be stated aS. (Fairbaim 1954), 

/■ (hV cos jgj) S dX = / (hV cos (p ) TdX + / (^- Dhv cos (p dX 

§ $ C 


+ / (i-~ DhU dcp- / hZ2^ds-iaa // fZcoscpd^PdX 

C L S 

+ ^ ;/ (FU + GV) cos<Pd(pdX (1.1,16) 

g s 

Note : We are carrying the frictional terms as far as we 
can even though our intention is not to retain them. This is 
done for the sake of a general formulation, and we drop their 
role at a convenient stage. 

In its ( tidal theorem) more practical form, if a certain 
family of free waves can be defined in the ocean area ACB 
but unrestricted by the boundaries, then the elevations and 
currents of forced tides along the parallel of latitude f 
can be determined by a series of explicit integrals involving 
the family of free waves, the tide generating potential, and 
the tidal elevations iround the coast ACB (Cartwright 1977). 

1 • 3 ' Auxiliary Functions 

For the sake of convenience we subdivide this section 
into two parts* The first part involves the auxiliary 



functions which are independent of the longitude coordinate x 


We call these auxiliary functions as auxiliary functions I, 
The second part involves the auxiliary functions containing 


the factor e 


ipX 


where p = 


n 




X and X ^ being 


the bounding longitudes of the region under consideration* 

We call them as auxiliary functions II, 

Auxiliary Functions I 

In this case the equations (1,1*8), (1,1*9) and (1*1,10) 
are written as 
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(1*3*1) 


Elimination of U leads to the system 
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where and Z^ are dimensionless and are given by 
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U, 
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(1.3.3) 


can be dete:mined from the relation 
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a cos cp 


(1.3.4) 


A uxiliary Functions II (containing a factor e^^^) 


Taking UyV^Z in the form given by 
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Z = aZ e^P^ 
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(1.3.5) 


where U / V and Z are dimensionless and functions of <P , 
P P P 

we can write the equations (1.1 .8), (1,1.9) and (I*! .10) in 


the form 
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Elimination of leads to the system 
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(1.3.8) 


1 Integral Formulation 


In the case of auxiliary functions I the integral 
formulation of the tidal theorem (li2.14) becomes 
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— zL— := o T^y 9 ^ y 

where ~ 2 ^ ^ H{e ^ — e and y.^ being 

the bounding latitudes for the region. Since the initial 
values for and are at our disposal, we choose 

($ ) = Z^($) =0 and V^(§) =1, Ihe last of the 
conditions is used in deriving equation (1.4.1). 

Similarly, for the auxiliary functions II with index -p, 
the integral fomulation of the tidal theorem (1.1 .14) becomes 
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Likewise for the index -p we have 
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where^ 

. - i'Tii __ i(p+2)X« i(p+2)Xi 

•^2 = 5 ^ ® ^ > 
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2-p 


- |Hi _ i(2-p)X_ 

e ^ H {e 
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and C = c + c . 
p p -p 

As the initial conditions are at our disposal we choose 

U, ^ (§) = O = Z (f) ^ V, (§) = i as initial values. The 
+ P iP ZP 2 

last condition is used in deriving the equation (1,4.4). 

Note ; In the equation (1.4,4) the surface integral containing 
Z is reduced to an ordinary integral/ by carrying out the 
integration with respect to X • C^/ c^/ c__p are the terms 

arisen out of the correction term S hZV ds. 

Method of determination of t along each successive parallels n 

Regarding the determination of the elevation ? along 
the various parallels of the region under study, the compu- 
tational procedure can be carried out in two different ways. 

In the first case, for the region ACB (Fig, Kb))/ we compute 
t along A3/ knowing the value along the boundary ACB, and then 
repeat the procedure for the region DABE to determine t along 
DE, Thus we proceed further downwards (Cartwright 1977). 



In this method, however, there is a continuous and e 3 <plicit 
dependence of g of the other parallels on S of the 
initial parallel AB (i.e,3 errors are carried). 

In the second case, we choose the lower most parallel 
of the region ACB (Fig, 1(c)) and determine t along AB, 
knowing the boundary values along ACB, We repeat this 
procedure for the region DCE to determine t along DE, 

Thus we exhaust all the parallels upwards. But this 
procedure suffers from the disadvantage of large magnitudes 
of the auxiliary functions due to large meridional extent of 
the region. We have chosen the second method for our compu- 
tations taking various difficulties of both the methods into 
consideration, 

Elevatiofi due to the constituent 

Let the extent of the bounding parallel § be from O 
to b. All the integrations in the previous section are 

■nTr 

carried out from 0 to b. Let p = ^ v^ere n is an integer 

Then the integrals are transformed to give the coefficients 

of the Fourier cosine series for t along , As n is ac 

our disposal, we can calculate the Fourier coefficients 

t 1 t f ^ etc. Thus with the transformation 
^o-' ^p' '*2p' 

b(X-X. ) 

X/ - - — we have . 

^ 2*^1 



59 


Sq - n ^ ^ ^^2,0 (l.S.l) 

and 5 = 2£ cos px' dx^ = S. + iS, „ (1.S.2) 

where refers to the real part and refers to the 

imaginary part. Having determined the mean elevation 
and the coefficients ^ 2p'' ^3p amplitude and 

the phase of the tidal constituent M 2 can be determined 
from the equation (1,1,4). That is from 


^M,= p’<'+2i,2p°°= 2px'+5j^_ 2pCOS 3pX' + ,*.)cos at 


•"it 2, 3pX''+..,)sin at 

(1,«.3) 


Tidal Data 


Tidal data# in the form of harmonic constants at 
various coastal stations are taken from the list of harmonic 
constants supplied by the Geodetic and Research Branch# 
Survey of India Department, The phase angles are given in 
terms of g 2 (given below) and necessary conversions are 
made to adjust them for a standarxi meridian and time. This 
is done using the formulae (Dronkers 1964) 


g2 = K + 2L - aT 
and gj^-g 2 s^ cr ( T-Tj^) 


(1 ,6,4) 



where < is the phase lag of the tidal constituent behind 
the phase of the corresponding equilibri^lm constituent at the 
place; L is the longitude of the station in degrees west 
of Greenwich; 

O is the angular speed of the constituent in degrees per 

mean solar hour (a = 28.984 deg/mean solar hour)y the time 
measured in standard time T hours later than Greenwhich Mean 
Time. 

The second of the above two equations (1.S.4) allows 
for a change in the time meridian frcm T to T^. Standard 
time is taken with T^ = -5.5 hours. That iS/ 5.5 hours fast 
with respect to Greenwich Mean Time* The phase lag is 

calculated for each station and these values together with 
the amplitudes are given in the following table. 

Table 1.1 


SI .No. 

Station 

phase iag 
in 

degrees 

Ampl itude 
in 

meters 

1 

Nagapatnam 

253.3 

0.1996 

2 

Madras 

237.3 

0.3316 

3 

Kakinada 

247.6 

0.4680 

4 

V is alchapatnam 

2 39.0 

0.4752 

5 

Paradeep 

241.4 

0.6218 

6 

Akyab 

254.0 

0.7821 

7 

Diamond Island 

254.42 

0^6712 

8 

Port Blair 

252.8 

0.6148 
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Use of this standard meridian, and the omission of the time 
factor gives. 



As the formalities regarding the formulation are almost 
complete, we proceed to describe the computational aspects 
in the following paragraphs. 


Co mput at ion of Auxiliary Function s 

For solving the systems (1.3.2) and (1.3,7) which are 

representative of the auxiliary functions I and II respectively 

a Runge-Kutta method based on Vemer's fifth and sixth order 

formulae has been used. In the case of the system (1.3 •2) 

the initial conditions at the bounding parallel § are 

specified as Uj^(f) = Z^(§) =0, V^(§) - ly and those for the 

system (1,3.7) are specified as U (§)=Z (f)=0, V = 0,5, 

jIP zp ■ zP 

The distribution of depth as a function of (P is fitted with 
the quadratic 

h ( cp ) = ( + a 2 ^ ) x 1 000 ,0 

where a = 1,360099, a- = 24,28704 and a^ = -72,42843, 
o 1 2 

The data for the bathymetry is taken from the charts provided 
by' the Naval Hydrographic Office » 



The extent of the latitude from the bounding parallel 

to the final parallel (i,e,# the open boundary at 20®N) has 

1 ® 

been partitioned into equally spaced strips of width j . 
Solution is advanced through each of these strips (in the 
direction of increasing ) using the Runge-Kutta-Vemer' s 
method* Solutions for 5 = 10,25 ®N are shown in 
Pigs, 3/4, 5/6 and 7* 

Interpolation of Harmonic Constants Along the Coast 

To evaluate the integrals on the right hand side of the 
integral formulation (1,1,14) we must know the values of Z 
along the coast-line which is available in the form of 
harmonic constants inferred from tidal observations at 
selected stations along the coast. It is necessary to 
interpolate than along the entire coast. For this purpose 
we have selected the rectangles of the mesh through which the 
coast— line passes (Fig, 2), There are 121 such rectangles 
(of size each) and the data is available for 8 rectangles 
only (Table 1,1), In such a situation as this, it is 
difficult to obtain proper distribution of the interpolated 
harmonic constants along the coast-line. From the available 
data for the harmonic constants for stations along the coast, 
it is observed that, except for the region at the head Bay 
of Bengal, the variations in phase lag are relatively small, 
but there is a continuous variation in amplitude fixxn 20 cms 
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near Nagapatnam to 75 cms near Akyab and slowly decreasing 
to 60 cms near Port Blair. The difficulty is further 
aggravated by the choice of an open boundary along the 
parallel at 20®N^ which has considerable extent^ at which 
the auxiliary functions attain their maximum magnitude and 
the amplitude of the tide too is quite high. The obtained 
tide is thus sensitive to the changes in the interpolated 
harmonic constants. Initially we have taken the linearly 
interpolated values of the harmonic constants {Table 1*2) in the 
rectangles along the coast between the coastal stations. 

However^ it became necessary to adjust the linearly interpo- 
lated constants to obtain a proper distribution of the tidal 
constituent along the parallel The adjustment is done by 

means of numerical experimentation - for the computation 
of the tidal distribution along the first parallel (Table 1*3) 
This adjustment is supposed to be valid for the computations 
regarding the tidal distribution for the other parallels as 
well. However^ we had to adjust the interpolated values of 
the harmonic constants once again for the reasons which will 
be mentioned in section 1.5, 

Correction Terms of the Open Boundaries 

If an ocean region has openings to the adjacent basins, 
one must account for the transport acres s these open boundaries. 
In the integral formulation (1,2.14) the transpoirt is 




Table 1,2 


Linearly Interpolated Values of Harmonic Constants 
G1 = PHASE-LAG IN DEGREES. H = AMPLITUDE IN METERS 


s. 

No. G1 

H 

1 

258.63 

0.156 

2 

256.85 

0.170 

3 

255 .08 

0.185 

4 

253.30 

0.200 

5 

251 ,5 2 

0.214 

6 

249 ,74 

0.229 

7 

247.97 

0.244 

8 

246.19 

0.25S 

9 

244.41 

0.273 

10 

242.63 

0.288 

11 

240.86 

0.302 

12 

2 39.08 

0.317 

13 

237.30 

0.332 

14 

237.30 

0.332 

15 

237.79 

0.338 

16 

238.28 

0 • 345 

17 

238.77 

0.351 

18 

239.26 

0.358 

19 

2 39 .75 

0.364 

20 

240.24 

0.371 

21 

240.78 

0.377 

22 

241.22 

0.384 

23 

241.71 

0.390 

24 

242.20 

0.397 

25 

242.70 

0.403 

26 

243.19 

0.409 

27 

243.68 

0.416 

28 

244.17 

0.422 

29 

244,66 

0.429 

30 

245 .45 

0.435 

31 

245.64 

0.442 

32 

246.13 

0.448 

33 

246.62 

0.455 

34 

247.11 

0.462 

35 

247.60 

0.468 

36 

247.60 

0.468 

37 

245.88 

0.469 

38 

244.16 ■ 

0,471 

39 

242.44 

0.472 

40 

240.72 

0.474 


S.No, G1 H 


41 

2 39 . 00 

0.475 

42 

2 39.00 

0,475 

43 

239.18 

0.486 

44 

239.37 

0.498 

45 

239.55 

0.509 

46 

239.74 

0.520 

47 

239 .92 

0.532 

48 

240.11 

0.543 

49 

240.29 

0.554 

50 

240.48 

0.565 

51 

240.66 

0.577 

52 

240.85 

0.588 

53 

241 .03 

0.599 

54 

241 ,22 

0.611 

55 

241.40 

0.622 

56 

241.40 

0.622 

57 

241.75 

0.627 

58 

242.10 

0.633 

59 

242.45 

0.639 

60 

242.80 

0.645 

61 

243.15 

0.650 

62 

243.50 

0.656 

63 

243.85 

0.662 

64 

244.20 

0.668 

65 

244.55 

0.673 

66 

244.90 

0,679 

67 

245.25 

0,685 

68 

245.60 

0.690 

69 

245.95 

0.696 

70 

246.30 

0.702 

71 

246.65 

0.708 

72 

247.00 

0.713 

73 

247.35 

0.719 

74 

247.70 

0.725 

75 

248.05 

0.731 

76 

248.40 

0.736 

77 

248.75 

0.742 

78 

249.10 

0.748 

79 

249.45 

0.753 

80 

249.80 

0.759 


S.No, 

G1 

H 

81 

250.15 

0.765 

82 

250.50 

0.771 

83 

260.85 

0.776 

84 

251.20 

0.782 

85 

251.20 

0.782 

86 

251 .20 

0.776 

87 

251.20 

0.769 

88 

251 .20 

0.763 

89 

251.20 

0.756 

90 

251.20 

0.749 

91 

251.20 

0.743 

92 

251 .20 

0.736 

93 

251.20 

0.730 

94 

251 .20 

0.723 

95 

251.20 

0.717 

96 

251.20 

0.710 

' 97 

251 .20 

0.704 

98 

251.20 

0.697 

99 

251.20 

0.691 

100 

251 .20 

0.684 

101 

251 .20 

0.678 

102 

251.20 

0.671 

103 

251.20 

0.671 

104 

251.29 

0.668 

105 

251.38 

0.665 

106 

251.47 

0.662 

107 

251.55 

0.65^! 

108 

251.64 

0.656 

109 

251.73 

0,652 

110 

251.82 

0.649 

111 

251.91 

0.646 

112 

252.00 

0.643 

113 

252.09 

0.640 

114 

25 2.18 

0.637 

115 

252.27 

0.634 

116 

252.89 

0.612 

117 

252.98 

0.608 

118 

253.07 

0.605 

119 

2.5 3 , 1 5 

0.602 

120 

253.24 

0.599 
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Table 1,3 

Adustment I of the intejrpolated harmonic constants 
G1 = PIIASE-LAG IN DEGREES. H = AMPLITUDE IN METERS 


s. 

No, 

G 1 

H 

s. 

No. 

G 1 

H 

s. 

No. 

G 1 

H 

1 

253.20 

0.195 

41 

239,70 

0.481 

81 

252.20 

0.735 

2 

252.30 

0.196 

42 

2 39.80 

0.486 

82 

252.70 

0.730 

3 

252.40 

0.197 

43 

240.00 

0.492 

83 

252.40 

0.725 

4 

252.00 

0.200 

44 

240.30 

0,500 

84 

252.00 

0.720 

5 

251.50 

0.215 

45 

240.60 

0.508 

85 

251.70 

0.718 

6 

249 ,50 

0.224 

45 

241.00 

0,516 

86 

251.40 ■ 

0.715 

7 

247.50 

0.234 

47 

241.50 

0.528 

87 

251.50 

0.712 

3 

245.50 

0.248 

48 

242.20 

0.540 

88 

251.80 

0.710 

9 

243.50 

0 . 2,62 

49 

243.00 

0.550 

89 

252.00 

0.708 

10 

242.50 

0.280 

50 

243.80 

0.560 

90 

252.00 

0.705 

11 

240.50 

0.295 

51 

244.70 

0.570 

91 

252.00 

0.703 

12 

238.50 

0.310 

52 

245.10 

0.576 

. 92 

252,00 

0.700 

13 

236.50 

0.332 

53 

245.10 

0.581 

93 

252 .20 

0.697 

14 

235 .50 

0.335 

54 

245.10 

0.537 

94 

252.40 

0.695 

15 

236,10 

0.338 

55 

245.10 

0.592 

95 

252.20 

0.693 

16 

236.50 

0.342 

5 6 

245.20 

0.595 

96 

252.20 

0.690 

17 

237.00 

0.346 

57 

245.20 

0.598 

97 

252 .00 

0.688 

18 

237.50 

0.350 

58 

245.20 

0.601 

98 

251.60 

0.685 

19 

238.00 

0.354 

59 

245.30 

0.605 

99 

251 .20 

0.632 

20 

2 38.40 

0.358 

60 

245.30 

0.610 

100 

251,00 

0.680 

21 

238.50 

0.362 

61 

245.40 

0.615 

101 

250.50 

0.678 

2 2 

238.80 

0.366 

62 

245.40 

0,619 

102 . 

250.00 

0.675 

23 

239.10 

0.370 

63 

245.40 

0.622 

103 

. 2 49.50 

0.670 

24 

240,20 

0,375 

64 

245.40 

0.625 

104 

249.9 0 

0.665 

25 

240.50 

0.385 

65 

245.50 

0.630 

105 

250.10 

0 . 660 

26 

241 .30 

0.395 

66 

245.50 

0.635 

106 

250..20 

0.655 

27 

241.70 

0.405 

67 

245 .50 

0.640 

107 

250.30 

0.650 

28 

241.90 

0.415 

68 

245.50 

0.645 

108 

250.40 

0.645 

29 

242.30 

0.425 

59 

245.60 

0.650 

109 

250.50 

0.640 

30 

242.70 

0.435 

70 

245.70 

0.655 

110 

250.60 

0.635 

31 

243,10 

0.435 

71 

246.00 

0.660 

111 • 

250.60 

0.630 

32 

243.80 

0.440 

72 

246.90 

0.665 

112 

250.70 

0,626 

33 

244.50 

0.450 

73 

247.70 

0.670 

113 

250.70 

0 ,62 3 

34 

245.00 

0.468 

74 

248,50 

0.680 

114 

250.80 

0.620 

35 

246.00 

0.460 

75 

249.20 

0,690 

115 

250.80 

0.618 

36 

247.10 

0.464 

76 

249.50 

0.700 

116 

250.90 

0.616 

37 

244,40 

0.468 ■ 

77 

249.90 

0.710 

117 

250.90 

0.615 

38 

242.10 

0.472 

78 

250.9 0 

0.720 

118 

251.00 

0.615 

39 

240.50 

0.475 

79 

251.00 

0.730 

119 

251.00 

0.613 

40 

2 38.10 

0.478 

80 

251.50 

0.740 

120 

251.40 

0.611 



represented by the integral / hZ2^ ds vAiich has to be 

OB 

evaluated along a fictitious line constituting the open 
boundary. For the evaluation of this Integral one should 
]cnow apriori the transport component v across the open 
boundary. For a general ocean it is difficult to obtain 
data regarding this, and there is no adequate data for the 
region under consideration. In the absence of such 
information an alternative representation of this integral 
is suggested by Fairbairh (1954). That is, when certain 
mean values for $ ,Z and <p are assumed viiich are represen- 
tative of the boundary enclosing the region S, then 


and hence 


/ hJ^ds= |r SSl dS 
OB 

/ hZ ds = Zp A 


(1.S.5) 


where A is the area of the basin adjacent to the open 
boundary. 

i(at - gg) 

with Sg = Hg e . we have 

A = a? cos // d<}> dX 

® A 

■■■ 


and 


( 1 . 5 . 7 ) 
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Since we have incorporated two open boundaries/ one 
towards north along the parallel at 20®N and the other along 
a fictitious line drawn connecting the submerged ridges of 
the Andaman Islands and the Diamond Island off Burma coast* 
Hence the correction terms in the integral formulations 
(1*4*1) and (1.4*4) will have two components which atrB 
representative of these open boundaries* ThuS/ in the 
integral formulation (1,4*1) we have 


C 


N 


n 

% ® 



- 

Hg e z 


E 


(1.5*8) 


Similarly for the integral formulation (1*4*4) with index p 
we have 




(1.5.9) 


1,6 Computations and Results 

A criterion has been specified in v;hich the matching 
of the computed value and the observed value of at 

the extremities (coastal points) of the bounding parallel 
is required. After the first adjustment of the linearly 



interpolated harmonic constants/ numerical experiments 

have been performed with the help of the above specified 

criterion* The matching of the computed values and the 

observed values for is achieved by adjusting the 

2 

correction tejan » This is because the only course 

left in the formulation/ after the adjustment of the 
harmonic constants/ is the adjustment of the correction 
terms to a proper level. The adjustment of p is 
done by choosing a proper combination of the mean values 
Hb, Zg/ Xg/ gg representative of the open boundary towards 
the north. It is the adjustment of Xg/ in most cases, 
which helped in achieving the required matching# The 
values of Hg and gg are hept the same as those for tiie 
case of computation for the first bounding parallel and 
the value of Xg is chosen as 20°N/ at which the values 
of Z_j_p are taken for computing the integrals. However, 
the value of Xgj^ has been adjusted in the range between 
91 °E and 93°E to achieve the required matching. The values 
of the integrals an<ifthe coefficients of the Fourier cosine 
series (in the equation (l,S,3)) are given in Tables 1*4 and 
1#5 respectively for the bounding parallel $ of 10®N* 



Table 


Integral 
/(S- Dv^dX 
-i /(S- 
/ IT dX 

G^J Z^cos^(pd(p 

S 

t = -0,1065 +0. 
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1,4 

Value 

0,013 + 0,003 i 

-0.005 + 0.096 i 

-0.053 + 0,037 i 

0.018 + 0,013 i 

0,002 * 0.006 i 

-0,002 - 0,013 i 

J34i 
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Table 1 .6 

Adjustment II of the Inte^rpolated Harmonic Constants 

Foir $ s 16 * 1 ^ ■ 

G1 = Phase-lag in Degrees. H = Amplitude in Meters 


S.NO. 

G1 

H 

S.NO, 

G1 

H 

S.NO. 

G1 

H 

1 

241.7 

0.405 

24 

244.4 

0.560 

47 

248.5 

0.675 

2 

241.9 

0.415 

25 

244.7 

0.570 

48 

249.0 

0,680 

3 

242.5 

0.428 

26 

245.1 

0.576 

49 

249,5 

0.695 

4 

242 .7 

0.438 

27 

245.7 

0.581 

50 

249 ,9 

0.710 

5 

243.1 

0.440 

28 

246.1 

0.587 

51 

250.9 

0.714 

6 

243.8 

0.440 

29 

245.1 

0,592 

52 

251.0 

0.720 

7 

244.5 

0.450 

• 30 

245.2- 

0.595 

53 

251.5 

0,728 

8 

245 .0 

0.468 

31 

245 .2 

0.598 

54 

252.2 

0.730 

9 

246.0 

0.473 

32 

245*2 

0.601 

55 

252.7 

0.730 

10 

247.1 

0.484 

33 

245.3 

0.605 

56 

252.4 

0.725 

11 

244.4 

0.490 

34 

245.3 

0,610 

57 

252.0 

0,720 

12- 

242.1 

0.502 

35 

245.4 

0.619 

58 

251.7 

0,718 

13 

240.-5 

0,500 

36 

245 .4 

0.622 

59 

251.4 

0.715 

14 

238.1 

0.510 

37 

245 .4 

0.625 

60 

251.5 

0.712 

15 

2 39,7 

0.514 

38 

245.5 

0.630 

61 

251.8 

0.710 

16 

239.8 

0.520 

39 

245,5 

0.635 

62 

252.0 

0.708 

17 

240.0 

0.522 

40 

245 .5 

0.640 

63 

252.0 

0.705 

18 

■240.3 

0.525 

41 

245 .5 

0.645 

64 

252,0 

0.703 

19 

240.6 

0.528 

42 

245.6 

0.650 

65 

252.0 

0.700 

20 

241 .0 

0.5 30 

43 

245.7 

0.655 

66 

252.2 

0,697 

21 

241.5 

0.5 35 

44 

246.0 

0.660 

67 

252.4 

0.695 

22 

242.2 

0.540 

45 

246.9 

0,665 

68 

252.2' 

0,693 

23 

243.0 

0,550 

46 

247.7 

0.670 

69 

252.2 

0.690 


% 


Matching of the computed tide and the observed tide 
(harmonic constants) at the extremities of the parallel § 
became a difficult task for higher parallels, particulatly 
after 16®N, At this stagej^ a second adjustment of the 
linearly interpolated harmonic constants became necessary* 
However, both these adjustments have been made within e 
tolerable difference. At any coastal rectangle the maxin'um 
phase angle difference due to adjustment with respect to the 
corresponding interpolated value is 10® , The corresponding 
maximum difference in the Case of amplitude is 5 cms# 
readjusted values of the harmonic constants are given in 
Table 1,6, This readjustment of the initially interpolated 
values of the harmonic constants has produced a reasonable 
matching of the amplitudes of the obtained and observed tides 
bringing down the highly, excessive amplitudes to an excess of 
10 cms above the observed amplitudes. In the following 
table we provide the values = 1,9, 3,4) for 

the latitudes 10®N, 12.5®N, 15?N, 16.25®N respectively 


Table 1.7 


Parallel 


10®N 


-0,01,0.27 0,27,-0.16 0.02,0.06 - 0 . 04 ,- 0.03 -O .01# -0-01 


12,5®N -0.15,0.31 0.23,-0.12 -0.06,0,03 0.04,0,05 




15°N -0,21,0.39 0.05,-0.18 0.03,0.08 -0.02#0,03— 0-01, -0.01 

16.25®N~0.22,0,43 0.08,-0.18 0.04,0.09 -0.03,-0.02 -O.Ol, 0,02 


73 


The longitudinal amplitude distributions for the 
latitudes 10°N, 12,5°N^ 15°N and 16.25°N are shown in Fig. 8. 

A cotidal map based on the tidal distribution on each 
of the parallels^ spaced apart, from lO^N to ie®N, is 
drawn (Fig. 9). 

The following concluding retffarks are made regarding 
the obtained results. 

Conclus ions « 

(i) Towards higher latitudes (above 18 ®n) the computed 
amplitude is 10 cm more than the observed amplitude. 

( ii) The computed phase lag is advanced by 45° with respect, 
to the observed phase. This phase lag remained 
almost the same in both the cases of adjustments of 
the harmonic constants about the initial interpolated 
values . 

( iii) to the case of readjusted values of the harmonic 

constants the amphidromic point (not seen in Fig. 9) 
is found shifted towards the west. 

( iv) Except for small discripancies in magnitudes the 

results have shown the trend similar to that found 
in the work of Pekeris and Accad (1969), 



List of Symbols and their Meaning 


6 

a the radius of the earth = 6. -37 x 10 meters 

to the angular speed of earth's rotation = 15 deg/hour 

g the acceleration due to gravity = 9*8 x 3600 x 3600 meters/hr 

h(^) the depth of water below any point in the mean surface 
tp the colatitude 

X the east longitude 

t the time measured from the instant of high-water of the 
equilibrium constituent at the meridian 80° E where the 
parallel of the latitude 10°N intersects the Indian 
coast of the Bay of Bengal 

2 (X/(P/t) the elevation of the free surface of the water at 
any time above any point of the mean surface 
S(Xrf<P^t) the equilibrium elevation corresponding to the t 

' i 

astronomical disturbing forces ; 

a the angular frequency of the tidal constituent* For 

M , a = 28*984 degrees per mean solar hour 
u,v,v the mean values along any vertical^ of the velocity J 

components at any time the directions of increasing X 
and ip and along the outward normal to a section 

[■ 

■I 

respectively 

i 

2 mathematical functions satisfying differential equations 
similar to those for tidal motion 
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N 


% 

H 


F,G 


the component along the outward normal to a 
boundary section of the vector whose components in 
the directions of increasing X and (p are U 
and V respectively 

the phase-lag at the place behind the phase of the 
corresponding equilibrium tide at Greenwhich 
the amplitude of the equilibrium harmonic consti- 
tuent, For H = 0,267 meters 

the components of exteacnal frictional force per 
horizontal area and per density of water 


P= 


7T 

X2~Xi 


5 



where Xj^ X 2 bounding longitudes 

enclosing the Bay of Bengal with the open parallel 
under consideration 

is the bounding parallel line in the open ocean 
along which we determine the distribution of the 
tidal constituent 

is taken for convenience in notation 
system of units is chosen for the entire 
investigation. 


M.K.H 
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FIG„2. MODEL AREA OF BAY OF BENGAL. 

i THE COASTAL STATIONS FOR WHICH THE 
HARMONIC CONSTANTS ARE AVAILABLE 
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CHAPTER 2 


DEPTH- INTEGRATED MODEL OF TIDES AND STORM SURGES 
IN BAY OP BENGAL 

2.1 Introduc-tlon 

Depth- integrated# two-dimensional models of ocean 
hydrodynamics are now very well established^ in which the 
basic hydrodynamical equations for the two-dimensional 
tides and surges are obtained by integrating the equations 
which are in full three-dimensional form/ from the bottom 
z = -h(x/y) to the free surface z = 5Cx,y/t). These 
models are called vertically integrated models. If the 
integrated equations are divided by the total depth (from 
z = -h(x/y) to z = $(x/y/t)) the models are known as 
vertically averaged models. These models, in the former 
or the latter form, have proved to be very useful in 
understanding the surge dynamics, particularly, for the 
prediction of sea-levels due to storms and other disturbances 
These models are relatively simple in conc^tion and require 
much less computer budget as compared to the esqjensive 
three-dimensional models. The equations, when reduced to 
complete linear form eliminate the vertical structure 
(Johns 1981). Some Important vertically- Integrated models 
have been mentioned in section 3, chapter 0, 



Though much literature on the subject for the region 
of North-Sea and elsevjhefe has appeared^ there is not much 
literature regarding the study of tides and surges for the 
region of Bay of Bengal, However there are a few important 
contributions to the study of surges for this region following 
the articles by Prank and Hussain (1971), Das (1972) and 
Flierl and Robinson (1972). Das et al.(1974) have studied 
storm surges including the interaction with tides and made 
a first major contribution on the subject for the region of 
Bay of Bengal. They have used the linearized long-wave 
eguations, dropping the role of advection and using a linear 
law for bottom friction with an assumed bottom current 
structure by a steady Elanan Spiral. Surges have been computed 
for three storm tracks (straightly directed towards North-West, 
North and North-east), using a Predictor-Corrector method. 

Johns (1981) has given a formulation for depth-averaged model, 
and a multi-level model with d— coordinate approach. In both 
the cases the atmospheric pressure-gradient is considered 
unimportant and neglected. Based on his approach Johns and 
All (1981) have studied the storm surges for the region of 
head Bay of Bengal and Johns et al. (1981) have studied storm 
surges for the east— coast of India. Recently Das (1981) 
has studied the tonporal response of wind— stress and has 

I hown that in the initial stages of the surge the divergence 
f the wind stress is more Important than the curl of the 
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wind stress^ while in the subsequent stages the curl becomes 
more Important. This response is related to the storm speed 
(vector motion) and is shown to be inversely proportional to 
the Square-root of the basin depth. An application of this 
model has been given by Sinha (1981) for the gulf of Thailand, 
Due to some limitation or the other^ many of these models^ 
though made a good start, lack some details needed for the 
better understanding of the surge mechanism. These models 
can be Improved by incorporating all the needed physical 
parameters into the problem, as many of these models have 
ignored one physical parameter or the other with the assumption 
of their negligible role. However, in particular local 
situations these physical parameters may be important and can 
no longer be ignored. For instance, the nonlinear advection 
may be important in extremely shallow waters and the so called 
* Inverted Barometer Effect' may be important in the case of 
slow moving storms (Dunn and Miller 1964), Also, in the 
nature, many of the storms are not axi- symmetric and the 
assumption of an axi-symmetric storm may considerably alter 
the magnitude, time and location of the peak-surge. On the 
other hand, the interaction with the tide is very crucial 
for the accurate surge prediction, A good model of surge 
requires the incorporation of all these factors. In view of 
the above factors it is felt that a reinvestigation of the 
surge problem with many needed details is essential. On the 
other hand the incorporation of essential features like the 



asymmetric wind field, the real bathymetry, the gradual 
Initialization and contirol of the storm movement etc», 
enhances our understanding of the surge mechanism paving 
the way for an accurate surge prediction. 

In the present chapter we have used the hydrodynamical 
equations derived in section 0.3 and incorporate all the 
essential parameters required for the study. In section 2 
of this chapter we have described the numerical model of 
the region under consideration. In sections 3,4 and 5 we 
have studied the three aspects of the surge problem, namely, 
the surge without tidal interaction, establishing of the 
tidal regime and the simulation of tide-surge interaction. 
Conclusions based on the results obtained in sections 3,4,5 
have been given in section 6. 

2.2 Numerical Scheme, and Model Area 

In the present chapter we consider the numerical 
scheme described by Flatter and Heaps (1975) in which the 
nonlinear advection is retained and properly represented 
by a scheme given by Robert and Wiess (1966). In the 
following paragraphs we describe the discretization 
procedure for the numerical treatment of the hydrodynamical 
equations governing the two-dimensional model of tides and 
surges derived in section 0*3. 
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An array of grid points in the (x^y) plane is taken/ 
consisting of ^ -points, u-points and v-points* The points 
of each type form a rectangular network of 1 columns and 
m rows with a square mesh of size s« Rows run parallel to 
the X— axis and colxamns parallel to the y-axis. Point by 
point/ relative to the 2-points , the u-points are displaced 
a distance s/2 in the positive x-direction and v-points 
a distance s/2 in the negative y-direction (Fig, 1(a)), 

All these points form a staggered grid system in which the 
points of each type are nxambered 

i = 1(1) 1 , 1+1(1) 21 / 21+1(1)31 /...,(m-l)l+l(l) ml 

Counting in the positive 3o-direction along each row, moving 
from one row to the next in the negative y-direction. If we . 
are to compute the values of U/V or $ at the points other 
than their own type the following averaging is used 

u at vwpolnt i : u. = 

V at u-point i ! = C 

where values of u,v at their respective 

points i. The following notation is used for denoting the 
values of the variables or parameters at the grid points. 



At u-point i : u = Pg = 5 

at v-point i : V = a X^ = X J ; 

at S-point i : ^ = ^ 

where are the depths at 5-point i. As mentioned earlier 

a simple averaging is used for the total depth 5 +h at u 
and v~points. Thus 

at u-point i ; d^^ = (h^+ 5^+ ^^+ 1 + ^i+i V 2 / 
at vwpolnt 1 .. = (h^+ 

F,G denote the wind stress components at the sea surface 
in directions of increasing X/y respectively and X 
denotes the bottom stress parameterization as described in 
section 0*3* 

Treatment of Advective Terms 

As we have decided to retain the nonlinear advective 
terms^ a numerical treatment of these terms requires special 
attention. A number of conditionally stable, explicit 
schemes are given by Lax and Wendroff Cl960), Crowley (1970), 
Sielecki and Wurtele (1970) and others. In the present 
nximerical schome for the simulation of two— dimensional tides 
and surges, we make use of the numerical scheme given by 
Robert and Wiess (1966). In the following paragraphs we 
briefly describe the scheme. For this purpose we consider 
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the following expressions 


3t 3 X 3 y 


(2. 2 pi) 


and ^ + u ^ + V ^ 
at 3 X 3 7 


(2.2.2) 


By the use of the selected staggered grid (Fig, 1(a)) 
the above expressions can be discretized, for odd time steps 
with i increasing, as follows 


u . ( t+x)-u.( t) u.(t+T)-u. -(t+x) u. (t)-u. (t) 


j] 


1 r’^-? U. ( t+T) -U , ( t+T) 

.If x-»m i—m-fl / i—m i 

2 ^ 2 ^ s 


) 


V (t)+v. .(t) u.(t)-u.. (t) 

1 x-fl A X x4ra 

^ V 


) ] (2.2,3) 






(■ 


) 


T- \ J 


. V. ( t+x)-v. (ttx) V*(t)-V. (t) 

- / . \ r 1 ^ 1 — m i .1 14m 

+ V^(t) {— -g— — + S 




(2.2.4) 



Similarly for even time steps with i decreasing we have 


V^(t4-T)-Vi(t ) 1 r ,V^^^Ct+T;)-v^(t+T) 

T 2 L 2 


l4in^ ^ "" i+l ^ ^ ’ 1 ^ ^ 




■)] 




- V . ( t+x) -V . ( t+T) V. (t)-V. (t) 

1 <• X i4m , i~m 1 


i] 


(2.2.5) 


u .(t+x)-u,(t) u . . (t+x)-u . ( t+x) u.(t)-u. -(t) 

T"^ I ^ 3 ■ : !] + 

+ 1 |- ^i(t)+V^^i(t) U^(t+T)-U^^(t+T) ^ 


V. _(t)+v. _.i(t) u. (t)-u. (t) 

4. l."*ra i**rn-fl ^ x**m 1. ^ j 


( 2 . 2 . 6 ) 


For the even time step we calculate v first and then u. 
The over bar in the above expressions denotes the average 
values given by 



("i-l + 2 "i + 


(2.2.7) 


With the above representation of the nonlinear advective terms 
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we write the discretized form of hydrodynamic al equations 
of section 3-chapter 0^ from (0.3.26) to (0,3.28) with 
forward-time and central— space differencing, as follows 


S^( t+T)-Sj(t) . 

7 = - I {ca^(t)u^(t)-d^^j^(t)u^_^(t) 

v^(t)} (2.2«8) 

For odd time— steps with increasing i, the equations of 
motion are descretized as 


u . ( t+T)-u . ( t) u.(t) 

T ® “ “fs — ^ ^i^ ^ “^i-l ^ ^ ] 


’• h L^^i-m^t)+v.,^^^(t)}£u^_^(t+T)-u^(t+T)} 


+ {v^(t)+v^^^(t)}{u^(t)-u^^(t)}] 


k u^(t+T) ^ 2 ,^, .._ 2 ,^^ t 1/2 


+ f v^(t) [u^(t)+v^(t) ] 


P.(t) 

+ Pd^ - I [ ] 


^ - Pl(t+T) ] 


( 2 #2 mS ) 



V. ( t+T)-v. (t) 

JL 






v,(t) 


- -f;— [ ^i.n,< t+t) -Vj^( t+T) +v^( t)-Vj^^( t)] t+t) 


k v. (t+T:) ^ 1/2 G.(t) 

— Atr~ + 4 


peTTtT 




(2.2.10) 

and for even time steps, with decreasing i (taking v first) 


V.(t+T)-V. (t) . 

? “ - b [(Ui(t)-Hi^^(t))(Vj^^j^(t«)-V^(t+T:5) 


+ (u^_^j^(t) tl^^^^(t)) (v^(t)-v^__^(t)) ] 


!ii£ 

2 s 


[ V . ( t+i:)— V . _ ( t+t) +v . ( t ) -V . 
^ X i+in 1 — m 2 


-fu^( t) 


i 

^.(t+T) ^ „ 1/2 G.(t) 

:fe — [u2(t)+v^(t)] + 


k VjCt+x) 

e 




(2.2.11) 



95 


u^(t+T)-u^(t) 

T " 


U^(t) 

2s ~ [ ^i+i ^ •Hi£( t:)-u^__^ ( t)3 


L. 

4s 


[ (v^(t)+v^_j_j^(t)) (u^(t+T:)-u^_^(t+Ty) 



Finally for solving the problem numerically we take 
the equations (2. 2.8) to (2«2«12) with appropriate initial 
and boundary conditionSy which we specify at appropriate 
stages in the following sections. In the next section we 
briefly describe the numerical model of Bay of Bengal 
representing its boundaries and other details such as 
bathymetry etc. 



Model area of Bay of Bengal 


Usually one is interested to study the surges in the 
shallow water regions such as the continental-shelf^ estuaries 
and river outlets, v^ere they are predominant and devastating* 
However, in storm-surge modeling, it is proper to include a 
much larger area encompassing the whole extent of the cyclone, 
consisting of deeper seas out side the continental- shelf as 
well. With this in view two model areas, a large area 
( Fig .1(c)) and a shelf-area (Fig, 1(b)), have been considered. 
Three probe runs have been carried out for both the models, 
with the storm forcing (other parameters set as in Table 2#l) 
along the tracks I, II, III of both the models. These probe 
runs have indicated a maximum-excess of elevation of 0,5 
meters of the large area model over the corresponding values 
of the shelf area model at the same coastal location in the 
case of positive surge. On the other hand a storm initia- 
1 ization to mature stage with its centre fixed at a point 
on the open boundary of the shelf-area model provided a 
reasonable closeness of the elevations with those correspon- 
ding to the large— area model. Weighing this fact on one 
hand and the envisaged extensive investigation with the 
1 imited resources for the computation on the other, it has 
been decided to take up the shelf-area model for the further 
investigation. 
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Bathymetry of the Shelf-Area Model 

As a compromise between the real and the ideal bathy- 
metries we have chosen the smoothened values of the bathymetry 
at the elevation points. In this process the Ganges Canyon 
is fairly smoothed (Figs, 2). Ihe smallest depth 

near the coast is 10 meters and the greatest depth at the 
open boundary is 300 meters. 

Storm Tracks 

Storm tracks for the present investigation have been 
chosen to represent those of typical severe cyclonic storms 
which occur in the region of the Bay of Bengal, Three- t 3 ^es 
of tracks^ making different angles of attack with the coast, 
have been chosen. Tracks I and II have the normal and 
parallel incidence respectively with the coast line at the 
time erf landfall (Figs, 3(a), (b)). The landfall location 
for the tracks I and II is set near Chittagong, The landfall 
location for track III is taken right to the Hoogly estuary. 
Track III makes normal incidence with the coast at the time 
of land-fall (Fig. 3(c)), 

2,3 Surge Without the Tide 

In this section, we intend to perform a number of 
numerical experiments to study the several aspects regarding 
the effects of atmospheric forcing on the sea surface. Wind 
stress as well as the sea-level barometric pressure are 
considered for this purpose. 


Boundary Conditions 


We prescribe the following boundary conditions vdiich 
are to be associated with the system of equations (2 •2. *8) 
to (2«2*l2)« In the discretized form they are 

u^(t)=0/ for the east and west coastal boundary points - 

for all t > 0 

v^(t)=0, for the northern boundary points of the coast - 

for all t>0 (2.3.1) 

and v^(t) = j[( ( t) } +d^(t)u^(t) 

- d. -(t)u. -(t) + e. „(t)v. ^(t)] (2.3.2) 

i— 1 i-1 i— m i— m ■' 

for the open boundary towards the south - for all t > O 

This condition follows from the continuity equation 

( 2 . 2 * 8 ). 

Asymmetric Storm 

In many practical situations of surge prediction we 
require three factors^ namely,^ the magnitude/ place and time 
of surge-event. Most of the studies on the subject deal 
with ideal storms having surface wind field with circular 

'"i ' 

syrcTOStry# However^ in real cases many stoimis deviate 



from this ideal circular-symmetry affecting the above three 
factors* Taking this fact into account asymmetry in the 
surface wind field is introduced in the form of deviation 
from the circular symmetry by using the formula (Chang 
and Anthes 1978) 

V^SY "= ® 1 0 5. < 2H 

(2.3.3) 

VgY and represent the magnitudes of the wind 

speeds in symmetric and asymmetric storms respectively. 

The argument 0 is taken anti-cloc3cwise from the right-hand 
direction perpendicular to the direction of movement of the 
storm-centre. As the intensity of storm may change during 
the course of its motion, the value of oc can be specified 
from time to time during the course of storm movement. The 
above formula indicates a stronger right-hand side wind 
field and a weaker left-hand side wind field with respect 
to the direction of storm movement in comparison with the 
symmetric s torm from which the asymmetry is obtained. This 
asymmetry is taken from the numerical experiments along 
with the symmetric storm for providing a comparative study 
of the response of the bay due to these two types of wind 
fields. For the present study the value of a is taken equal 
to 0,3 for the entire course of storm movement. 



Time Step for Numerical Computations 


The governing inequality for the linear stability 


T < 


'/gh 


max 


(2*3.4) 


for the chosen grid space s ~ 13 km yielded a maximum value 
for time step i:=240 seo^nds. For overall stability in the 
computations a time step t ==120 seconds has been chosen 
after a probe.' run. 

Meteorological Forcing 


The following meteorological forcing parameters are 
considered for the numerical ej^periments 

Table 2.1 


Storm parameter 


Value 


Central-pressure drop; aP 


Radius of maximum winds 


mn 


Radial extent of the storm; R. 


max 


Radius of the eye ; R 


eye 


Speed of the maximiim sustained 

Average vector motion of the 

storm 

Angle of inflow in the annular 
region between R^^^ and 

at theimatured stage 


50 mill! bars = 5 x lO"^ Newtons/: 

sq. mtr, i 

52 kn 

I 

300 kn I 

I 

13 kn ^ 

60 meteirs/sec | 

i 

■ ■ ■ '■ i 

20 kr/hr 


20 ® 





Cyclostrophic wind profile 

formula 

Wind-pressure relation 

Drag coefficient for wind 

stress 


(Das et al#«1974) 

V = 0.85 '/’Ap 
max ^ 

= 0,003 (Simons 1974) 


Also we consider the following. 


Drag coefficient for bottom Ic - 0,0025 

stress 

Growth time to matured stage =10 hours 

of the storm 


Numerical Experiments 

We have performed two numerical experiments, with the 
help of the parameters described in the preceding paragraphs 
and table 2.1, In the following paragraphs we describe the 
numerical experiments and discuss the results obtained from 

these experiments. 

Experiment 1 . 

Experiment though an elementary one, is intended to get a 
first hand information regarding the steady— state response 
of the sea basin under study, due to an Imposed wind— stress 
field of a stationary storm whose centre is stationed at 
3^(20,19) on the open boundary of the model area, TOiis 
experiment is repeated by stationing the storm— centre at 
32 ( 20 , 32 ) of the model area (Fig, 1(b)). A point near 
Chittagong is chosen for noting the hourly elevations. Also, 



the hourly values of the maxlmuin elevation with respect to 
the entire shelf ( for each hour) have been noted* This has 
provided a valuable information regarding the shifting of 
the hourly maximum— surge in space for the region* Finally/ 
the elevations for the entire shelf area are recorded at 
the end of 35 hours of integration which was enough duration 
for establishing a steady-state response of the sea due to 
stationary storm* The e35>eriment has been performed with 
and without the role of advection* 

Results . 

This experiment has provided a testing case for the 
validity of the present model and an Insight into the order 
of magnitude of the elevations due to the prevailing meteoro— 
Ipgical forcing* The sea-level contours at the end of 
35 hours of numerical integration are shown in Figs* 4 (a) /(b)* 
The positive surge due to on-shore and the negative surge 
due to off-shore winds are clearly seen* This experiment has 
also shown the negligible role of advection/ owing to the 
stationary storm with centre situated far away from the 
coastal boundaries* In all the runS/ with storm centres on 
the open boundary/ the peak— surge elevations never exceeded 
1 meter* The magnitude of the storm forcing and the 
magnitudes of elevations are in reasonable agreoBent with 
the observations* This justifies the validity of the present 
model for carrying out further numerical experiments* 



Experimenta 


In the following numerical experiment^ unlike 
the previous experiment, we consider a moving storm along 
the tracks as shown in Figs. 3(a), (b), (c). Two sets of 
numerical experiments are chosen corresponding to the 
symmetric and asymmetric storms. Thus, keeping the other 
parameters fixed, we consider the following combination of 
parameters for the present numerical experiments which 
includes nonlinear advection. 

Symmetric storm ^^Track I 

along c;;^ — Track II 

Asymmetric storm-^ ^^Track III 

Probe runs have been carried out with the following 
storm durations. 

Period of runs for tracks I, II and III ; 60 hours 

Hour of landfall for tracks I and II : 35 hour 

Hour of landfall for track III : 26^^ hour 

During the first 10 hours of each run the storm is 
gradually allowed to reach the mature stage and then allowed 
to move along the chosen track. This has been done in order 
to avoid any spurious oscillations due to suddenly forced 
strong winds. For the first three runs of the e>^eriment a 
syiranetric storm moving along the tracks I, II and III 
respectively is considered. After the landfall the storm 



movement is continued further for a period of 10 hours 
during which time the storm intensity is brought down 
gradually to a 10 knot wind (1 knot = 0.5144 meters/sec.) . 
This procedure is adopted after carrying out the probe runs 
to determine the locations of maximum surge-peaks’'" for each 
of these runs in which no significant surge peaks have been 
obtained after the storm intensity has been decreased to a 
10 knot wind speed-after the landfall. The three subsequent 
runs have been repeated with the asymmetric storm as 
specified in the parameter list (table 2.1). 


The results computed for the runs with symmetric and 
asymmetric storms along tracks I, II, III include the hourly 
values of the elevation-maximum in space and its corresponding 
location. Also computed are the hourly values of elevation 
at the location of landfall, as well as at the location of 
another potential surge-peak, and the values of the individual 
terms in the hydrodynamical equations. 


The distribution of the sea-level elevations, due to 
asymmetric storm along track I, at the three selected 
instants of time (at 29, 35, 40 hours respectively) are • 
shown in Figs. 5(a), (b), (c); while those corresponding 
to the symmetric storm along track I are shown in 5(d),(e),(f) 
respectively. The corresponding elevation contours 


Maximum surge peak: Maximum taken 

for all time 
steps 


Maximum 
for all 
space points 




(at 29, 35/ 40 hours) for track II are shown in Figs* 6(a)/ 
(b)/ (c) and Figs. 6(d)/ (e), (f) respectively. The 
elevation contours due to asyrnmetric storm forcing along 
track III (at 21, 26, 31 hours respectively) are shown in 
Figs, 7(a)/ (b)/ (c) and those corresponding to symmetric 
storm forcing (along track. Ill) are shown in Figs. 7(d)/ 

( e) / (f) respectively. 

With asymmetric storm along track I the computed peak- 
surge of 4.8 meters, at the point of landfall, has occured 
3 hours before landfall while the computed maximum peak-surge 
of 8,6 meters has occured at the point (2,36), seven hours 
before landfall (Fig, 8(a)), The corresponding results 
with asymmetric storm along track II revealed a sustained 
negative surge at the point of landfall almost throughout 
the period of storm movement, and the maximum peak-surge 
of 6,9 meters at (13,46), ten'hours before landfall 
(Fig, 8(b)), Similarly the corresponding ccxnputed results 
for track III revealed a peak-surge of 2 ,0 meters at the 
point of landfall, ten hgurs before the landfall and the 
maxim\jm peak-surge of 3 ,5 meters, three hours before the 
landfall (Fig, 8(c)). These aspects are summarized in the 
following table. 



Table 2.2 


Asymmetric 
storm along 

Elevation above 
the Mean Sea 
Level in meters 

Location 

coordinates 

Time of occurence 
with respect to 
land fall time 

Track I 

MPS 

8*6 

(2,36) 

7 hours before the 
LFH 


SLF 

00 

• 

(3,43) 

3 hours before the 
LFH 

Track II 

MPS 

6.9 

(13,46) 

10 hours before the 
LFH 


SLF 

-3.0 

(3,43) 

8 hours before the 
LFH 

Track III 

MPS 

3.5 

(8,22) 

3 hours before the 
LFH 


SLF 

2.0 

(8,18) 

10 hours before the 
LFH 


MPSs Maximuin Peak— Surge? 

SLF: Peak— surge at landfall location? 

LFHs Landfall Hour. 


The numerical experiments with symmetric storm have 
shown that the peak-surge elevation can be 2 .0 meters less 
than those corresponding to those with the asymmetric storm 
(Figs.- 8(a), (d))^ However the corresponding maximum peak- 
surge locations, for any track, remained almost the same in 
both the numerical experiments* The computed results 
corresponding to the synmetric storm are sunmarized in the , 


following table* 




Table 2*.3 


Symmetric 
storm along 

Elevation above 
the Mean Sea 
Level in meters 

Location 

coordinates 

Time of occurence 
with respect to the 
landfall hour 


MPS 

6.7 

(2,36) 

8 hours 

before the 

Track I 

— 



LFH 



SLF 

3.5 

(3,43) 

3 hours 
LFH 

before the 

Track II 

MPS 

4.8 

(13,46) 

8 hours 
LFH 

before the 


SLF -3.0 

(3,43) 

7 hodrs 
LFH 

before the 

Track III 

MPS 

2*6 

(8,22) 

4 hours 
LFH 

before the 


SLF 

1.3 

(8,18) 

14 hours before the 
LFH 


where the abbreviations are as those given for Table 2 •2*. 


In all the above mentioned numerical experiments the 
maximum peak-surge as well as the peaky^surge at the landfall 
point have occured several hours before landfall* An interes- 
ting observation from the present numerical experiments is the 
combined role of the coastline configuration of the bay and 
the angle of attack of the storm near the coast, giving an 
interesting behaviour of the surge patterns* Ihe maximum 
peak-surge to the left of track I (normal to the coast) and 
the maximum peak— surge to the right of track II are the mani- 
festations of such a track-ciastline combination (Figs*. 8 (a) 
(b) and 8(d), (e))* 


9 




In the case of numerical ej^erlments with track III, 
however, there are some local, short-lived, abrupt oscilla- 
tions around the point of landfall (Figs. SCc) and 8(f)). 

This has resulted in sharp fall— rise— fall of th^' elevation 
with sharp peak at t = 22 hours near the location of land- 
fall. This is also evident from the depth— mean current 
vectors depicted in Figs. 8(g) to (1) from t = 12 hours to 
t = 22 hours at two— hour intervals. From these figures we 
obsejTve the anptying by the westward current and the subse- 
quent swift filling by Northward current at the point of 
land— fall. The large horizontal extent of the coast and the 
approach of the storm in the normal direction are presumed 
for this kind of swift oscillations. This is particularly 
so because of rapid change in the wind-direction near the 
point of land-fall as the storm closely approaches the coast. 
This kind of situation is peculiar to this kind of coast- 
track combination and is not observed in the case of numerical 
experiments with track I and track II. 

The individual contribution of the terms in the hydro- 
dynamical equations, in the form of vectors (also called stick 
diagrams) at hourly inteirvals, at the three selected locations 
for each of the tracks I, II, III (with asymmetric storm) 
are shown in Figs.. 9(a), (b), (c)? 10(a), (b), (c) ; and 
11(a), (b), (c) respectively. In the figures these are 
termed as CORIOLIS, i^VECTION, AIMOSPHERIC PRESSURE, WIND 
STRESS, ELEVATION SLOPE, BOTTCM FRICTION respectively. In 



each case the vectors are nonnallzed with respect to the 
maxinrum of the magnitudes of the vectors of all the terms 
taken together. From the results we observe that the 
contribution of the atmospheric pressure (the so called 
Inverted Barometer Effect) is insignificant in all the cases. 

We also observe that in most of the cases of peak surge^ the 
elevation slope has played a major role, particularly during 
the hours before landfall. It is the restoration force of 
the surface slope following the bottom contours, particularly 
in shallow areas, vhich caused the majority of these peak 
surges. The set up by the wind is rather a gradual accxmu- 
lative process (except when the storm centre approaching the 
landfall point) as compared to the local surface slope 
dynamics ( free- wave types) , However, during landfall hours 
the direct action of the atmospheric forcing may be in 
'resonance^ with the forces due to the surface slope 
inequalities. In such cases the surges can be more devastating 
The swift local oscillation in the case of track Ill- 
numerical experiment is, perhaps, due to such manifestation 
of resonance. However, in shallow areas such a resonance is 
controlled by the bottom friction, which has comparable 
magnitude to counter it - bringing such an oscillation 
confined to local scales in space and time. This can be 
seen in Fig . 10 (a). On the other hand, for the region under 
consideration, the nonlinear advective terms played no 
significant role in the surge dynamics - even in the shallowest 



areas upward. This envisages the inclusion of a fuirther 
shallow^river system in which nonlinear advection can be 
significant. In all the cases it is observed that the 
depth— mean currents as well as the elevations decrease 
rapidly with the decreasing intensity of the meteorological 
forcing after the landfall, 

2.4 Establishing the Tidal Regime 

As a first step for the study of tide-surge interaction 
we establish an appropriate tidal regime giving steady 
CO— oscillations of the tide corresponding to a prescribed 
sinusoidal oscillation of 1 meter amplitude, in phase, along 
the open boundary of the region under study. For this purpose 
we consider the equations (2 •2,8) - (2 *•2 #12 ), with the 
exception that we drop the meteorological forcing 
GI = O) and include the tidal forcing as a boundary condition 
across the open boundary. Thus the associated boundary 
conditions for the present problem, in the discretized form are 

= 0 for the east, west coastal boundary points for t > O# ^ 

0 for the northern coastal boundary points for t > 0, 

=! A cos 0 t , t > 0, where a denotes frequency of the 

M^-tide and A denotes the amplitude of the ^ 

total tide, ' 

.(2.4 a) 

for the open boundaiy points f®^ t ^ ^ 


u^( t) 
v^( t) 
S^Ct) 

v^Ct) 



The last equation follows from the continuity equation. 
This equation will make use of the elevations prescribed as 
tidal forcing at the elevation points adjacent and interior to 
the open boundary • In addition to the above conditions we 

have 

=Vj^= ^^=0 at t = 0 for all the interior points of 

the region* In the present case we have taken A = 1 at all 
the points of the open boundary. 

Two values for the time— step are chosen satisfying the 
necessary condition for the stability of the numerical scheme 


T < 




(2*4*2) 


max 


and probe runs are carried out with t = ^ Lunar hour and 


T = ^ Lunar hour. 


It was found T = ^2 hour (116 


seconds approximately) yielded satisfactory results. With this 
time step we have performed two runs corresponding to the 
hydrodynamic al equations with and without advective terms. A 
criterion is set for establishing a satisfactory tidal regime 
by which the integration process is terminated, if at the end 
of a cycle the elevation at any point is close to that with 
respect the corresponding value at the end of the previous 
cycle, within a tolerable difference of 1 cm* In the case of 
no advection the tidal regime with the desired accuracy is ' % 
achieved after 28 cycl^ of Integration* In the run with 
advection the desired tidal regime is achieved after 30 cycles 



of integration* However the tidal regimes achieved in hoth 
the cases are very close and nonlinear advection did not play 
an important role in the tidal dynamics for the region under 
study* 

The elevation contours are given in Figs* I2(a ) to 12( 1 ) 
and the d^th-mean currents are given in Figs* 13(a) to 13(1)* 

The computed tidal regime has produced a High Water of amplitude 
2*3 meters as compared to the observed amplitude of 2 meters 
near Chittagong* The computed phase lag of High Water near 
Chittagong with respect to the High Water near the open boundary 
is 120 degrees ( ~ 4 hours)* These computed values are 

fairly in agreement with the observed values* Depth-mean current 
ellipses for the eight selected locations are given in Figs* 15 (a) 
to (h) , At shallow area comers the depth-mean tidal currents 
have shown distorted patterns indicating their complexity at 
such locations (Fig* 15(3)), Ihis is also evident frcxn the 
current vector diagrams (figs* 13(e) and 13(j))* For a ccanpa- 
rison on time-scale of the elevation and the depth-mean current* 
tidal profiles are provided in Figs* 14(a), to (h) at the eight 
selected locations* The model has responded well to the 
presence of the 'swatch of no land' (Ganges Canyon) making the 
main crest of the tide to refract and progress in two directions^ 
North-West and Noridi— east (Figs* 12(a ) to 12( 1 ))• 



2 *5 Tide— Surge Interaction 


It is well known that surges and tides interact with each 
other modifying each other's propagation* In shallow areas 
like an estuary this interaction may be important in view of 
the nonlinear mass transport^ and for a good prediction of 
surge this factor must be taken into account* In recent years 
some progress has been made regarding the study of tide-surge 
interaction and a general account on the development of this 
subject is presented in section 0,2, For the region of the 
Bay of Bengal, tide-surge interaction has been studied by 
Das et al* (1974), Johns (1981) and Johns and Ali (1981). In 
the present investigation we have incorporated some of the 
features not considered by these earlier studies* For instance, 
we have included the nonlinear advection terms which are omitted 
by Das et al. (1974). The Importance of these terms can be 
seen in the case of surge due to meteorological forcing along 
track I (Pig* 9(a). In contrast to the works of Johns (1981) 
and Johns and Ali (1981) the present study differs in the 
selection of parameters such as storm tracks, wind-field 
asymmetry. 

Model . 

We have chosen the same model as has been described in 
the section 2*2 regarding the meteorological forrcing. An 
additional feature is the specification of tidal forcing- 



The main feature 


as an Initial-Boundary Condition"'^ ( IBC) , 
of the meteorological forcing is the asymmetric storm along 
tracks 1/ II and III, 

Tidal Forcing . 

We have selected the tidal regime established in the 
previous section for the purpose of astronomical tidal forcing. 


By Initial-Boundary condition we mean the specification 
of the values of u^v, ^ corresponding to a particular phase 
of the established tidal cycle over the entire region at 
time t = 0 (initional condition) / and the subsequent con- 
tinuation of the sinusoidal tidal oscillation along the open 
boundary for further time steps (Boundary condition). 

Thus 

u. = V. = V. tide, ?• = t = 0 for all i 

rail IX 

and = 1.0 X [cos (dt)]^ for t > 0 at the open boundary 

points 

and ~ "e ' "rt) ' ^ ^ ^j[^t)u^(t) 

for t > 0 at the open boundary points. 


( 3 . 5 . 1 ) 




For each track we have performed 4 numerical experiments with 
4 types of initial— boundary conditions representing the four 
phases of the established tidal cycle corresponding to the 
Maximum (High-Water of 1 meter) — Falling Zero — Minimum 
Rising Zero of the elevation along the open boundary*. 

As in the case of pure-surge prcblem (sec* 3*2) we have 

out the probe runs and found no significant displacement 
of the maximum surge— peak locations as those corresponding to 
the pure-surge problem* Taking this fact into account we have 
retained the same locations for noting the elevations as those 
corresponding to the pure-surge problem. The surge computations 
are carried out with the IBC referring to the phase values along 
the open boundary and the corresponding phase lag for any 
coastal location can be found from the figures 12(a ) to 12( !)*• 
All the numerical experiments are performed through 4 cycles of 
integration approximately and the obtained results are plotted 
in Figs, 16(a) to (d) , 17(a) to (d), 18(a) to (d) respectively. 

In contrast to the pure surge problem where, at any 
coastal location, we normally have single surge-peak, several 
surge peaks are obtained, each corresponding to the tidal cycle 
at the location. In almost all the cases the linear supers 
position of the tide on the mirge produced higher magnitudes 
of peak surge than the corresponding surge peaks due to tide- 
surge interaction. The exception is the phase 3— track II 



case where the peak surge at the location (8^22) coincided 
with the Lov^-Water phase at the location* However, at any 
instant of time 8 location, the difference between the elevation 
due to tide— surge Interaction and the elevation due to tide+ 
surge is relatively small ( fv o,5 meters) as compared to the 
magnitudes of the corresponding elevations thanselves* Also, 
the positive surge-peaks, at a location in most of the cases 
under investigation have coincided with the High Water at the 
location* These factors indicate that the tide-surge inte- 
raction is rather weak in the present model area and the 
inclusion of a further shallow area into the river system 
would have been more proper as the trends of the present results 
indicated a stronger interaction in the estuar^river system 
which will be a large enough extent for the propagation and 
interaction of tide and surge*. 

2*6 Concluding Remarks 

Based on the results of the nximerical experiments, we 
would like to spell out the following concluding remarks, 

l) The specification of the wind field is necessary 
rather than an idealization from the given pressure drop* 

Many of the cyclonic storms possess some kind of asyanmetry 
regarding pressure and wind field distribution* The 
specification of the pres sure- wind relations in terms of the 
radial distance r (fron the centre Of t^e storm) becomes 



difficult in such cases / and as a result we have many such 
relations in practice* We have shown that for a given 
radial pressure distribution of the symmetric storm the 
derived asymmetry produced more devastating surges-with 
the magnitudes differing by as much as 2 meters in some cases* 
With the present kind of monitoring by satellites and other 
devices we should be able to specify the surface winds more 
accurately. Efforts should be concentrated in this regard. 

2 ) Coastal geometry/ bottom contours and the storm 
tracks play a vital role in the surge dynamics. It has been 
found that the Bay of Bengal has all these features (peculiar 
to itself) whose combination in surge dynamics are quite 
interesting. For instance/ at the same land fall location/ 
tracks I and II have produced entirely different kinds of 
surge response with different locations/ time and magnitudes 
of maximum peak surge* On the other hand, track III has 
produced a different kind of surge response with local/ 
short-lived, swift oscillations near the point of land fall- 
during the land fall hours. All these are attributed to the 
track-coastline- bottom contour combination* In the case of 
the numerical experiment with track III, the obtained unusual 
behaviour is believed to have been caused by some kind of 
resonance due to local effects* On the other hand it is 
also believed that this oscillation is confined to the 



surface layer of that local region (owing the quick damping) 
Perhaps the current structure through the depth may provide 
a better answer for this problem. In the next chapter we 
explore this aspect. 
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CHAPTER 3 


NUMERICAL SIMULATION OF THREE-DIMENSIONAL 
SURGE USING A SPECTRAL MODEL 


^ J-ntrod uctlon 


In certain situations it will be useful and interesting 
to know the vertical structure of the ocean characteristics 
such as the ocean currents and the ocean themocline* Though 
the vertically-averaged models are quite powerful in the 
studies dealing with the elevations and the depth-mean 
currents, the three-dimensional models are more realistic 
and provide very useful information regarding the vertical 
structure of the ocean, particularly when the ocean circula- 
tion is predominant in the vertical. In sections 0.2 and 
0,4 we have mentioned some of the well known three-dimensional 
models. Many of these models are quite recent and there 


are relatively a few number of studies as compared to their 
two-dimensional counterpOUts , The reason is simple, of 
course. Many of the studies either do not require them or 
avoid them because of the e 3 q)ensive involvanent. But in 
the studies of residual circulations, and other circulations 
these models play a vital role. Also, it can be seen from 
the literature that the ipcsaei'^et .extensive studies, using 


the three-dimensional imodels, have begun only in the 1970s, 
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Many o.€ these models follow either the multilevel (or multi- 
layer) approach or the spectral approach. In the present 
work we consider the spectral model given by Heaps (1976) for 
the purpose of simulating the three-dimensional surge in the 
Bay of Bengal, This choice provides us some of the advantages 
of the spectral model over the multi-level model. The 
apparent advantages are, the easy and natural way of represen- 
ting the bottom topography; the direct way of satisfying the 
dynamical conditions at the free surface and at the bottom; 
the determination of the currents at any desired level (as 
opposed to the determination of the currents at fixed levels 
only) and the achieving of sufficient accuracy in the deter^ 
mination of the vertical structure by increasing the number 
of eigen functions, and hence the number of terms in the 
summations. We consider the same kind of meteorological 
forcing as has been taken in Chapter 2, regarding the tracks 
and other storm parameters. However, only asytanetric storm 
will be considered because of the limited number of numerical 
experiments that can be carried out. As in the previous 
chapter# the storm forcing will be t 3 ^ical of a severe cyclonic 
stom in the Bay of Bengal, 

# 

As regards the simulation of surges for the region of the 
Bay of Bengal, to the best of the author^ s knowledge, only one 
significant study has been carried out with the help of a 
multi-level model (Johns et al. 1932(a)), The high-light of 



the present investigation is the discussion on the vertical 
structure and the residual circulation patterns. However, 
a comparison of the results of the present work with the work 
of Johns et al. (1982) is not possible because of the 
entirely different features of the works. The only anchor 
is the set of few available observations for the region of 
the Bay of Bengal, where some kind of comparison can be made. 

In the following section we briefly describe the numerical 
scheme and then proceed to the aspects of numerical experiments. 

3*2 Numerical Scheme 

In the present investigation we adopt the numerical 
scheme given by Sielecki and Wurtele (1970) for solving the 
system of equations (0.4.31), (0.4.32) and (0,4.34) with the 
appropriate initial and boundary conditions* The same scheme 
has also been considered by Heaps (1976) in his formulation to 
include the nonlinear terms. With the use of the staggered 
grid taken in Chapter 2, we have (for the variables ^,u^,v^ 

9-t S-«point i i , h = hj_ 

at u-point i : n=u^, ^r" s i 

at v^point i : i* i'^s'^i 

simple averaging yields esttaates for total depth at u-point 1 
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( 3.2.2) 


Similarly the total depth at v-point i is given by 


e. = (h^ H- J. + h.^ + 


(3.2.3) 


In the same way other averaged eguantities are taken as 
follows; 

u at v-point i : u. = (u. .+u,, )/4 

^ 1 1-1 1 i+m-1 i-Hna ' 

V at u-point i : v^ = ^’^i-m’*'^i-m+l'^^i‘^’^i+l^^^ 

"r v-point i! 

Vr u-point i: v^,i=(Vr^l,V^r,i-n.+l+^r,r'^r,i+l'/'* - 


(3.2.4) 

Also# on the basis of the ofuadratic law for bottom friction 
we have 


. u k 
X • * 


‘r, i c 




V 


? t(S )f.+(v,)f^)^/Ve^ 


■r,l - S “r "“i'-h^'-i 


(3.2.5) 


where from (0.4.33) 


( 11 ) ss T ^ u S COS O* 

^^i -h *r rj-i ^ 

(\r ) as V ^ V ^ 

^^i -h r,i r 


(3.2.6) 



Introducing a time step t and including the terms 

^U/t * r ^v,r' proposed numerical scheme 

by Sielecki and VJurtele (1970) can be used to write equations 
(0.4. 34)/ (0.4.31) and (0,4.32) respectively in the finite 
difference form as follows: 


g^(t+T)-?^(t) 


M 

«r' -ai.i ( ( t) 


+e. (t)v . (t)-e.(t)v -(t)}/s 

i-m r, i-m 1 r/i 

(3.2.7) 


i^ t-f t)~u^^ j^( t) 




ga. 


{ ( t+T) - Si( t+T) ( t+T) +?<( t-f-T ) } 


F^(t) 


+ R . - S 

Uv.r, 1 u^r, X 


(3.2.8) 


v„ .(t+T)-v .(t) 

.X,#..;X. , JL, , 


a - .(t) V (t+T)-fu ^(t+ir) 

r j 1 T / 1 r. i 


ga 


£ c ? , ( t+T ) - t+-t)-q( t+T) +q^( t+T) > 


G . ( t) 


)e^(t) ^v,r,i + 


(3.2.9) 


followed by a second step of the form 
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'"r,! '^r,i<t+T;’ + f? i(t+T) 


ga. 


T“ ^ ^i4.l('t+T)-?j.(t+T)-^^_j_j^(t+T)+?^( t+T) } 


F.(t) 

+ \~k I \.\ + R 4 - S ' 

pa.(t) u,r>i u,r 4 i 


(3*2,10) 




= -’ .(t+T)-fu“ .(t+x) 


‘r^i^ r^i 


r^i 


ga 


^ { ?i(t+i:)^q^^(t+T)-q(t+T)-r._^(t+x)3 


Gi(t) 

"'' p^TTtT ^v,rf± “ ^v,r,i 


(3.2.11) 


where = ( V 1-l^r, l-*r,i-Kn-l+*^^, 


(3.2.12) 


The cjuantities ^^(t+x), u^^i^t+x), vj^^(t+x) represent 

the final values of the fields at t+x as deduced fTOm those 

at time t using (3.2,7) to (3.2.11) in succession. In the 

e< 5 uations (3.2.8) to (3.2.11), the finite-difference forms of the 

nonlin^sT* tsms R. #R ,S -S • S' and S_ are 
liuniinear rexms V 4 .r' n,r* v,r' u,r v, r 


n 




(3.2 .13) 



V, r, 1 e 


j_ ( t+x)-^ .( t) 

JtT ^ ^ ^ 


H-v .(t)+ S (1+a b .)f .V. .(t)} 

r r,j 3 3,1 


(3.2.14) 
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ji n 


+ r) v.(i,i+l,t) 6u (l,i4TO,t)] D 

J ^ ^ / J / Ia 


•S. ri 


V, r,l 
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{A,(t)+/i,^(t)! J (3.2.16) 
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4d*^Xt^ S 2 W j^! j ( i**! # i# t) j ( i# i+1 # t) } 


{Ai(t+-r)+A^^^(t+T)) ^ (3.2.17) 
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where 
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(3.2.19) 




and 



In all the alx)ve sunmations, the indices j and n are carried 
from 1 to M, The expressions for B . d , and E 

r/jn 

can be found from (0.4.65) to (0,4.69). 

Applied successively^ in the order stated/ these 
relations may be used explicitly to advance from the fields 
of i time t to the corresponding fields 

at time t+T / thereby forming the basis of a marching 
technique for computing the quantities ^ .^u . and v . 

1 V f 1 XT/ 1- 

through space and time. 

Boundary Conditions 

we 

In the present three-dimensional model, j_cons ider the 
surge without the tidal forcing. Owing to the limitation 
on the computer budget, the fresh establishment of the tidal 
regime (which is necessary) by means of the spectral model 
was not possible. Under this situation we specify the following 
appropriate boundary conditions. 

At the coastal boundaries we get 

u . tsa V . ts 0 for t ^ 0 (3.2.21) 

ITiif J. JT# J. ' . ' 

At the open boundary we set 

for t > 0, 


q = 0 


(3.2.2g) 
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v^v ,„v "f-'T C 


h .1 W 

"i l,i l,i 


n 


+ x ih%“ .'jI ,(t)-h'; ,a" . 

1 r# i T/ i r# 1 i-l r / i-l r, ,i— 1 r i— 1 


t > 0 (3.2.2 3) 


and . =0, t > 0 , ' 

lE # .1 


(3.2 .24) 


We set the initial valties of and equal to 

zero, which correspond to the undisturbed sea. With the 
above initial and boundary conditions and with the specified 
storm forxiing, we seek the response of the sea region through 
the numerical solutions. Knowing the values of i ^ ^r, i 
and the current structure is determined by the relations 
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u, « E « S cos(a §) 
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(3,2.25) 


The depth-averaged currents are given by 
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(3.2.26) 


currents near the 




-)‘v: 



fr'i-e- surf ace, 
givc-'D by 



and 



Thus, the currents at the free surface are 



n u . 
r r.L 


(3.2 .27) 
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Similarly, setting = -1 in (3.2.25) the bottom currents 
(which are reqfuired for the bottom stress parameterization) 
can be obtained. 


3.3 Nu merical Experiments 

In the present investigation, the model shelf area, 
the asymmetric storm forcing and the tracks (Fig. 3(a), 
(b), (c) of Chapter 2) are kept the same as in the case of 
tvro-dimensional surge simulation (Chapter 2). The main 
features of the storm forcing arc given in Table 2.1 
Surge-tide interaction is kept out of the scope of the 
present chapter because the model requires the generation 
of a fresh tidal regime which should provide all the 
spectral modes at any required phase of the tidal cycle. 
The establishing of a fresh tidal regime requires a large 
computer budget and constrained by the limited resources 
the study is confined to the steuiatio« of surge alone. As 
in the previous chapter a pjx^* carried out to 



extract first hand information regarding various aspects 
such as the location# time and magnitude of the maximum 
pGaT< surciGj the number of terms (modes) required for a 
reasonable convergence of solutions etc,- The grid spacing 
is taken as in the case of two-dimensional surge (Chapter 2 ) 
and the corresponding time step of 2 minutes is taken to 
satisfy the overall stability criterion (2.3.4), The 
probe run is carried out with the asymmetric storm forcing 

along track I for a period 46 hours taking the landfall 

til 

at 35 hour. With the results of Chapter 2 in view# 
regarding the weak role of the nonlinear advection in the 
model# the advection terms are included. The value of ^ 
in the equation (0,4.54) is taken equal to o■o^o^^ , , 

which represents an average value for the present model. 

The first five modes are taken as significant modes and the 
rcanaining modes are neglected. 

The results of the probe run have indicated the 
behaviour# regarding elevations# similar to that of the 
corresponding two-dimensional model. However# the 
elevations obtained in the spectral model are slightly 
higher in magnitude as compared to the elevations of the 
corresponding run of the two-dimensional model. The 
overall weak role of the nonlinear advection is also seen 
in the spectral model, . ; as well as the prohi 

bitively huge computer budi^ of the model has 



coiTipcllC'cl the modol to drop the advection terms in the 

subsequent mns. 

Witli the above considerations we have min the present 
rnocic'l three times^ with the storm forcinq alonq the three 
chosen tracks. Based on the information from the probe 
runs the following plan has been envisaged for recording 
the results from the model executions, 

1) For each run# corresponding to the storm forcing 
along the chosen track# throe locations# namely# the 
location of maximum peak surge# the location of land fall 
and the location of a second maximum peak surge (same as 
those of the corresponding run of the two-dimensional model) 
have been selected for noting the hourly elevations. The 
coordinates of these locations for the runs for each track 
are given below 

Table 3,1 

Coordinates of the selected locations 



Location of Location of 
MPS LF 


Location of 
second MPS 



9) Ten locations have been chosen to examine the 
current, structure through the depth. At these locations 


the v-i.1ufS of u,v and the roots of the equation (0.4.54) 
at hourly intervals are recorded for the first five eigen 
values (corresponding to the first five modes), with 
these values the current structure at any required level 
in the range between ^ = -1 (bottom) and g = 0 (free 

surface) can be determined using the equations (0,4.33) . 
The selecte<1 locations include the locations mentioned in 
Table J.l. Ihese locations are shown in Fig, 18 whose 
coordinates are given below. 


A1 (?,i6) 
A7 (4,40) 
Ai (6,28) 
A4 (3,16) 
A5 (8,20) 


A6 (10, 16) 
A7 (11, 20) 
AS (13, 24) 
A9 (13, 46) 
A10(18, 6 ) 


For the pueposB of presenting the current vectors and 
the elevation pattern for the entire shelf area under 
consideration, three time intervals have been chosen. These 
results, for each run, are presented at three instants of 
time which are given below 


Surface curren 


Bottom current 


Depth-moan current 



Four hours before landfall 


At the time of landfall 
"Npive hours after landfall 


1 IrMl ly the hourly maximum peak surge is noted to have 
a comoAratlvc look with respect to the corresponding hourly 
iar.:diriur; (in space) peak surge of the two-dimensional model. 

'^?ith this plan of noting the values^ the runs of the 
model arc executed with the duration of the storm forcing 
along each track as taken in the corresponding two-dimensional 
model. 'Itsat is# 



track 

I s 

46 

hours 

with 

the 

landfall 

at 

2gth 

hour 

along 

track 

II s 

46 

hours 

with 

the 

landfall 

at 

3^th 

hour 


track 

III s 


hours 

with 

the 

landfall 

at 

26^^ 

hour. 


Results and Discussion 

All the runs of the spectral model have revealed slightly 
large magnitudes of the surge elevations as compared to the 
corresponding surge elevations of their two-dimensional 
counter parts. This is clearly shown in Figs. 1(a) through 
Kc), 2(a) through 2(c) and 3(a) through 3(c). in most 
of tlio cases any significant difference in magnitudes of the 
corresponding elevations at any instant never exceeded 1 meter 
Also, tho corresponding evolution patterns of the elevation 
at a location remained almost the same in both the models. 
Figures 4(a) through 4(c), 5(a) through 5(c) and 6(a) 
through 6(c) show the elevatixsn : dontour bands at the three 



selected instants of 



storm torclnq along each track). Except for small difference 
in magnltiJ'icSif these contours are similar to the corresponding 
con ton rn f)r the two-dlnensional model. As in the case of 
two-dimensional model for track II, the landfall location in 
the precent model has experienced negative surge around the 
hours ol landfall. Fig. 7 shows a comparative behaviour of 
the elevation pattern due to storm forcing along tracks I 
and II at the point of landfall. Owing to the orientation 
of track II, the storm has, most of the time during its 
movoment, produced off-shore winds at the point of landfall. 

On the other hand, the storm along track I has produced 
north-ward long-shore winds causing a positive surge at the 
same point of landfall. The evolution of the elevation 
patterns at the three selected locations (for each track) 
are shown in Figs. 8(a) through 8(c). These figures can be 
compared with the Figs. 8(a), (b), (c) of Chapter 2. Except 
for the excessive magnitudes of the elevations near the 
concave comer points, all the above results indicate a 
close similarity, to a large extent, in the elevation patterns 
of the dopth-averagod model and the spectral model. 


Unlike the elevation patterns, which present smoothly 
varying gross features, the ocean currents are difficult to 
conceive and analyze. This is particularly so, far the 
surface currents which gulekly fesipnd to the strong winds 
and for which no reliable 



current mctc'r rc-idinqs very near to the surface are not 
avail aM<, evn-rs for the moderate wind conditions. Also^ the 
attrlHition of tJic current structure is difficult (in the 
absence; of observations) because the resulting current 
donondc not only on the direct external forcing but also 
on the ovi.'rall dyiKimics. Taking these limitations into 
account, the residual current patterns are presented for 
the region under consideration under the prevailing strong 
wind forcing. 


Ibc current vectors at the free surface, at the bottom, 
and the depth-mean current vectors for the numerical experiment 
with the storm forcing along track I, at the three chosen 
instants of time, arc shown in Figs. 9(a) through 9(c)i 
10(a) through 10(c)j and 11(a) through 11(c) respectively. 

The corresponding current vectors for the numerical experiment 
with the storm forcing along track II are shown in Pigs. 12(a) 
enrough 12(c) j 13(a) through 13(c) J and 14(a) through 14(c) 

rospoctivoly. similarly the corresponding current vectors 

for tho numerical experiment with the storm forcing along 
track III arc shown in Figs* 15(a) through I5(c)j 16(a) 
through 16(c)| and 17(a) through 17(c) respectively. 


In all these figures we observe the horizontal circulation 

wi-ch an inflow at the central .{ortion of the open boundary, 

, 1-1 Ai^^if^fw^tlal-shelf (shown as 

and out flow side ways alo?^ v 



'IN' and 'O-JT' in tho figures), Ihis is expected because 

oJ tilt.- -/idc extent of the basin (comparable to the size 

of the storm), and the smaller depth of the basin as 

compared to the so called 'Ekman depth' (E.E = Ji ^2h/I ) 

A is the coefficient 
of vertical friction 

for th<- region under consideration ( s ee Appendix) . This is 
found to b<' in agreement with the horizontal circulation 
pattern described by Welander (1961). Perhaps the comparable 
sizes of the bay and the storm (Harris 1957, 1959) might 
have resulted in this type of circulation. We observe the 
concentration of outgoing surface currents across the ends 
of the open boundary. In all the numerical experiments 
the strongest currents prevailed to the south-west of the 
storm centre and this strong residual circulation prevailed 
many hours after the storm crossed the coast. 


Another major feature is the large magnitudes of the 
surface currents (more than 6 meters/sec. \ i ) in the 
core wind region during the development stage of the 
stationcxl cyclone (Figs, 22(c), (d) and 24(c), (d)). 
However, when the storm has started its movement these 


surface currents have shown a decreasing trend in magnitude. 
Long fetch of accelerating surface winds upto mature stage 


of tho cyclone has produced stronger surface currents than 
the maximum winds of the cyclone. However, the 
persisting strong surface-' shelf en 's o 



the 01HJ3 to\indary are believed due to the strengthening of 
the r-H rent:; by the winds and bottom slope. Figures 25, 

9 6, 97 ( '1 t.o <'1) Show that the bottom currents to a large 
extent have the s<iine direction as that of the surface 
currentr. chowinq either (locally) weak or no vertical 
ciroilnLion in the region of study. However, at the present 
no reliable field observations are available regarding the 
wind driven currents for the region of the Bay of Bengal 
for a comparison of the present values. If one disregards 
the period of the developing stage of the cyclone, the 
present values (magnitudes) of the surface currents, when 
conpared with the studies with similar storm forcing, are 
in agreeable range. 


Figure 31(d) shows an abrupt oscillation of the 
v-component of the current in the surface layer during the 
landfall hours while there is a continuous change of the 
magnitude of the u»ccxnponent fron negative to positive 
during the sot© period. This fact and the presence of the 
coast towards the north are believed to be the cause for 
the swift accumulation of the water mass in the vicinity 


(particularly to the right) of the point of landfall. 
Howevor, as anticipated in chapter^, this oscillation is 
confined to the surface layer (top |th). Except for this 
swift adjustment the overall variation in the elevation 
(or current) is smooth and kind of 



o.sciHaHon is a very particular situation. In fact it 
has b'.cn verified with the help of the two-dimensional 
model/ that no swift oscillation has occured if we choose 
the track obliquely in either direction (passing through 
th<' same landfall point of track III). 

Conc lud inq Remark s 

I'ho numerical experiments in the present chapter have 
dutionstratcd the variety in the bay response to the storm 
forcing. From the llmiited number of results taken from the 
three numerical ej^erljnents (which has resulted in a 
paltry 61 figures i ) the following conclusions are made, 

1) Regartilng elevations/ both the models (Depth- 
averaged and spectral) have produced similar kind of 
response. However/ the elevations obtained in the three- 
dimensional model have shown higher magnitudes (up to 10 % ) 
than the. corresponding elevations of the two-dimensional 
model. It is assumed that the formulation of the bottom 
stress with the bottom current and the limited number of 
terms (modes) taken in the three-dimensional model might 
have resulted in this discripancy* 

2) Under the strong storm forcing the bay has responded 
with a predominant horizontal circulation with an inflow 
near the central portion- tii@ <m 1 ^ 1 ,Qw sideways across 


the open boundary. 



i) rh<’ model has shown strong surface current in the 
surliic'o 1 I'/' T during the developing period of the cyclone. 

At Gornr local ions these currents have shown magnitudes- 
in i;xc(;r5 of the normally expected values, A varification 
of this was not possible at the present stage^ due to 
non-availability of observational data, 

4) The local# swift oscillation of the surface elevation 
duo to the storm forcing along track III is explained in a 
better way as comoared to the explanation given in the depth 
averac-H'id model# showing the worthiness of the three- 
dljnens Jonal model in particular situations. It is also 
shown that the oscillation is a very particular case and is 
confined to the surface layer (resulting in quick damping). , 
No swift oscillation has been observed in the numerical 
experiments with the storm forcing taken along the oblique 
tracks (either way) passing tdirough the same point of 


landfall. 
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' CHAPTER 4 


A CRITICAL REVIEW AND SOME SUGGESTIONS 


We conclude the present dissertation with a critical 
review of the work carried out in Chapters 1^2 and 3, Some 
suggestions will be made regarding the improvanent of the 
present work and the inclusion of new aspects such as the 
coastal upwelling and the thcimnal structure of the Bay of 
Bengal as a response to the translating cyclones and the 
large scale wind circulations such as the monsoons* 

In the present dissertation an attempt has been made 
to study some details of the dynamics of tides and surges 
in the Bay of Bengal thaxiugh numerical experiments* The 
study has been carried out in two stages. In the first 
stage (Chapter 1) , the distribution of M^-tide is determined 
for the region of the Bay of Bengal (from 10°N to 20®N) 
using the tidal theorem enunciated by Proudman (1925), The 
availability of the limited niomber of coastal- station 
observations (harmonic constants) made the present study to 
rely much on numerical experimentation through matching of 
the obtained distribution and the harmonic constants at the 
coastal stations* The recent acquisition of some additioiteJ. 
tide guages ( to be set up along the coast) would help 



further investigation regarding the determination of the 
open-sea tides for the region of the Bay of Bengal, The 
co-tidal map obtained in the present work (Chapter 1), for 
the region of the Bay of Bengal, is reasonably close to the 
observations and agrees with the trends indicated by the 
co-tidal chart in the work of Pekeris and Accad (1969), 

In the second stage (Chapters 2 and 3) surges, tides 
and their interaction are studied with the help of a 
vertically averaged model (Chapter 2) given by Flather and 
Heaps (1975), and a spectral model regarding the three- 
dimensional surge (Chapter 3) given by Heaps (197G), The 
moving coastal boundaries are not incorporated in the 
formulation of the surge problem, and as a result the 
magnitudes of the elevations (and currents) are found to be 
in excess of the observed values. At the present stage 
a comparison of the results of the present work, regarding 
elevations and currents, with the field observations could 
not be made because no adequate and reliable field data is 
available for the region of the Bay of Bengal. The recent 
acquisition of the research vessel ('Sagar Kanya' by the 
National Institute of Oceanography at Goa, India) , the 
geosynchronus satellite (INSAT IB) and the setting up of 
additional tide guages along the coast etc., would, perhaps, 
provide a base for the field data for the future investigations 
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Owing to the nature of the shelf area at the head 
Bay of Bengal, many complex types of motion are possible. 

Some of these are the swift, local oscillations obtained 
in the numerical experiments with track III (Chapters 2 and 3) 
and the entirely different behaviour of the surge at the 
same point of landfall due to storms along different tracks. 
The models used in Chapters 2 and 3 have shown a relatively 
weak response regarding the nonlinear advection and the 
interaction between the tide and the surge. These nonlinear 
physical processes are believed to be stronger at smaller 
scales (of space and time) and hence require a greater 
resolution regarding the grid spacing, particularly in the 
shallow areas. The inclusion of further shallow areas 
such as the river mouths and the nearby Islands in the model 
area is necessary for the actual surge prediction. 

The tidal regime established in Chapter 2 has produced 
tidal oscillations which are remarkably close (regarding 
the amplitudes ) to the observed tides near the coastal 
stations. But the depth-mean current (as indicated by the 
current ellipses) has shown some complexity with the rotation 
of the current vector in the clockwise direction at some 
locations. At the present a justification of this is not 
possible due to lack of adeguate field observations for the 
region under study. This task is left for a future 
investigation. 



The depth- averaged model as well as the spectral model 
of the present work have shown remarkable similarity regarding 
the elevation patterns and the depth-mean current structure. 

An interesting result is the horizontal circulation under 
the strong meteorological forcing where the storm has the 
horizontal extent comparable with that of the shelf at the 
head Bay of Bengal. The present investigation has considered 
only the strong meteorological forcing typical of severe 
Cyclomic storms in the Bay of Bengal, It will be interesting 
to study the circulation patterns under moderate, large-scale 
wind conditions such as the monsoon winds. The baroclinic 
response and the coastal upwelling have to be taken up in the 
future investigations extending the scope and purview of the 
present models. 

The present study also points out the need for adequate 
field observations for a comparison of the theoretical and 
obseirvational studies. For instance, the numerical experi- 
ments of Chapter 3 indicated a strong surface currents during 
the developing (intensification) stage of the cyclonic 
storm ( even during movement) and a varification of this 
aspect is not possible owing to the nonavailability of the 
pertinent information regarding the surface currents. 

Excepting the surface currents at some locations the overall 
current structure through the depth are v/ithin the physically 
possible (or observed) range. 
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With the results of the present work in view^ the 
following problems are suggested for future investigations: 

i) the role of the large submarine canyon (the 
Ganges Canyon) in the dynamics of tides and surges in the 
Bay of Bengal, 

ii) the barotropic and the baroclinic response of 
the Bay of Bengal to the translating cyclones* In the 
recent studies it is observed that the development and the 
movement of cyclones are sensitive to the changes in the 

sea surface temperatures. As a part of air-sea Interaction, 
it will be interesting to study this aspect in which the 
thermal structure of the ocean and the storm (and its 
movement) interact with each other. This study helps the 
prediction of the course of the storm movement and hence 
the prediction of the surge. With this view, consideration 
should be given to the large-area model (which we have 
discarded in Chapter 2) in which the entire course of the 
storm movement can be incorporated* 

iii) circulation over the shelf and the coastal 
upwelling due to cyclonic storms and large scale wind 
conditions such as the monsoons. 

In spite of the inclusion of the aspects such as 
the real bathymetry, the asynmetric storm forcing, the 



nonlinear advection# the tide-surge interaction and the 
three-dimensional current structure, the present work has 
achieved the objectives in a ' limited way» What has been 
achieved in the present investigation is short of a tiny 
droplet. The droplet is showing the abundance of the 


ocean 



APPBITDIX 


•RT rman Deipth for the Region of Bay of Bengal 

Wind effects penetrate to some considerable depth 
below the ocean surface, determined in part by the stability 
of the water column. The usual limit for wind effects (such 
as drift currents) is taken to be the depth at which the 
current arrow is exactly opposite to the surface-current 
arrow. At that depth, the current speed is about 4 percent 
of this surface current. Under a strong wind, wind- drift 
currents have been observed to extend to depths of about 
100 meters. The depth of the Ekman Layer is given by 

E.D. = n Jzk/f, ■ 

where A is the coefficient of vertical friction and f is 
the coriolis parameter. In the present investigation we 
take the value of f = 0-5 x \o^ which corresponds to the 
latitude 20° E and A^q.qts wS- an average value for the 
Bay of Bengal. Thus the Ekman depth for the region under 
consideration is meters. Also we that in most part 

of the region under study we have the depth of the basis 
less than this ’ Ekman Depth’* 
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